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ABSTRACT 

We present a gravitational wave (GW) analysis of an extensive series of three-dimensional magnetohydrodynamical core-collapse 
simulations. Our 25 models are based on a 15M Q progenitor stemming from (i) stellar evolution calculations, (ii) a spherically sym- 
metric effective general relativistic potential, either the Lattimer-Swesty (with three possible compressibilities) or the Shen equation 
of state for hot, dense matter, and (iii) a neutrino parametrisation scheme that is accurate until about 5ms postbounce. For three 
representative models, we also included long-term neutrino physics by means of a leakage scheme, which is based on partial im- 
plementation of the isotropic diffusion source approximation (IDS A). We systematically investigated the effects of the equation of 
state, the initial rotation rate, and both the toroidal and the poloidal magnetic fields on the GW signature. We stress the importance 
of including of postbounce neutrino physics, since it quantitatively alters the GW signature. Slowly rotating models, or those that do 
not rotate at all, show GW emission caused by prompt and proto-neutron star (PNS) convection. Moreover, the signal stemming from 
prompt convection allows for the distinction between the two different nuclear equations of state indirectly by different properties of 
the fluid instabilities. For simulations with moderate or even fast rotation rates, we only find the axisymmetric type I wave signature 
at core bounce. In line with recent results, we could confirm that the maximum GW amplitude scales roughly linearly with the ratio of 
i-C ■ rotational to gravitational energy at core bounce below a threshold value of about 10%. We point out that models set up with an initial 

central angular velocity of 27r rads 1 or faster show nonaxisymmetric narrow-band GW radiation during the postbounce phase. This 
emission process is caused by a low T/\W\ dynamical instability. Apart from these two points, we show that it is generally very diffi- 
cult to discern the effects of the individual features of the input physics in a GW signal from a rotating core-collapse supernova that 
can be attributed unambiguously to a specific model. Weak magnetic fields do not notably influence the dynamical evolution of the 
core and thus the GW emission. However, for strong initial poloidal magnetic fields (> 10 12 G), the combined action of flux-freezing 
and field winding leads to conditions where the ratio of magnetic field pressure to matter pressure reaches about unity which leads to 
the onset of a jet-like supernova explosion. The collimated bipolar out-stream of matter is then reflected in the emission of a type IV 
GW signal. In contradiction to axisymmetric simulations, we find evidence that nonaxisymmetric fluid modes can counteract or even 
suppress jet formation for models with strong initial toroidal magnetic fields. The results of models with continued neutrino emission 
show that including of the deleptonisation during the postbounce phase is an indispensable issue for the quantitative prediction of 
GWs from core-collapse supernovae, because it can alter the GW amplitude up to a factor of 10 compared to a pure hydrodynamical 
■ treatment. Our collapse simulations indicate that corresponding events in our Galaxy would be detectable either by LIGO, if the 

source is rotating, or at least by the advanced LIGO detector, if it is not or only slowly rotating. 
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1 . Introduction lescence from black-holes (BH) and neutron stars up to distances 

of several 100 megaparsecs, but might also very likely produce 

Gravitational wave (GW) astronomy may soon become a reality the first detections of black-hole and neutron star mergers, core- 

and will allow humankind to address questions about many dif- co ii apse supernovae, a nd neutron star normal mode osci llations, 

, ferent astrophysical objects that are hidden from the electromag- as rece ntly reviewed in Sathvapra kash & Schutzl I2TJ091) . 

■ netic detection. Within the past few years, the first generation of xhe broad ma j ority of GW sources can be subdivided into 

the ground-based GW detectors LIGO (USA), VIRGO (Italy), two classes: (i) mat hematically well-posed events such as BH- 

GEO600 (Germany), and TAMA (Japan) have got very close to BH coalesC ence and (ii) scenarios that involve matter. The first 

or even reached design sensitivity and collected partially coinci- categor y of sources can be modelled accurately, the waveform 

dent data, as discussed by| Abbott et al. | (|2009j). Lately, the two 4 can be calculated with high precision, and thus a matched-filter 

km LIGO detectors were upgraded to sensitivities increased by analysis can be applied . Qn me other han(1 the latter kind of 

a factor of 2-3 (enhanced LIGO, see |Adhikari et al.| (2006)) and GW sources C an only be modelled imperfectly, as matter effects 

resumed observations in 2009. Upgrades of the three LIGO in- carry large physics uncer tainties and open up a huge parame- 

terferometers and VIRGO are expected to be completed by 2014 ter space for initial conditions . Moreover, even if there were 

and will increase the observable volume by a factor of ~ 1000. < per fect' models, there is still a turbulent, stochastic element 

Operating at such a high level of precision, GW detectors are in some GW emission mechanisms that makes it strictly im- 

sensitive to many different sources, such as compact binary coa- possib i e t0 compute templates. As a result, more general, un- 

modelled bu r st ana lysis techniques are used, as summarised in 

Send offprint requests to: S. Scheidegger lAbbott et alJ {2009)' and references therein. One of the scenar- 
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ios where burst data analysis must be applied is the stellar core 
collapse, where in addition to the matter effects even the funda- 
mental explosion mechanism is not fully settled. It is still un- 
clear which processes may convert some of the released grav- 
itational binding energy of order (9(10 53 )erg into the typically 
observed kinetic and internal energy of the ejecta of (9(10 51 )erg 
= 1 Bethe [B] needed to blow off the stellar envelope. While 
most state-of-the art simulations investigate a mechanism based 
on neutrino hea ting in combina t ion wi th hydrodynamical insta 
bilities, such as Ma rek & Jankal d2009l), others explore the alter- 
native magneto-rota t ional mechanism dLeblanc & Wils'onl 



1970; 



2008 



Kotake et al. 2004a; Burrows et ID 120071 Mikami etal 
Taki waki et al.l2009l) . or t he acoustic mechanism whi ch is driven 
by strong core g -modes (Burrows et al.ll2006l 120071). For a re - 
cent review see Uanka et all (120071) iTakahara & Satol (fl 988). 
Gentile et all dl993h . and lately also iNakazato et alJ (120081) and 
Sagert et all (|2009) reported that a QCD phase transition may 



power a secondary shock wave which triggers a successful hy- 
drodynamical explosion. 

Beside neutrinos, which have already b een observed in th e 
context of stellar core collapse of SN1987A dHirata et al.l fl988>. 
GWs could provide access to the electromagnetically hidden 
compact inner core of some such cataclysmic events. Since a 
core-collapse supern ova is expected to show aspherical features 
dLeonard et ali 2006). there is reasonable hope that a tiny amount 
of the released binding energy will be emitted as GWs which 
could then provide us with v aluable information about the angu 
lar momentum distribution ( Dimmelmeier et al 



baryonic equation of state (EoS) jMarek et a" 



2008 



200' 



and the 
both of 

which are uncertain. Furthermore, they might help to constrain 
theoretically predicted SN mechanisms (Ott 2009b). Given the 
actual se nsitivities of today 's operating ground-based GW obser- 
vatories (IWhitcombll2008h combined with knowledge of wave- 
forms from state-of-the art modelling, GWs from a Galactic 
core-collapse supernova may be considered to be within the 
detector limits. The collapsing iron core and the subsequently 
newly formed PNS can be subject to a whole variety of asym- 
metric hydrodynamical and nonaxisymmetric instabilities which 
give rise to GW emission. 

In this context probably most attention during the past three 
decades has been paid to rotational core collapse and bounce 
dynamics. In this particular phase, the core spins up due to 
angular momentum conservation, resulting in an oblate and 
time-dependent deformation that leads t o strong GW emission 
(Mull erl 119821: iMqnchmever et alJll99lb iJanka & Mullerl[ l996: 



Zwerger & Mullerlll997URampp et alJI 19981: iKotake et alj| 2003. 
2004bi lO bergaulin geretalJ 120061: lOttetalJ 120041 l2007allbfc 



Dimmelm eier et al.ll2002ll2007Ll2008tT Scheideg ger et al.l l2008). 



The steady improvement of the models (e.g. the inclusion of GR, 
a micro-physical EoS, treatment of neutrino physics) recently 
led to a theoretically well understood single and generic, so- 
called type I wave form which is characterised by a large nega- 
tive peak at cor e bounce, f ollowed b y ring-down oscillations that 
damp quickly dOtt et all d2007allbl) : IDimmelmeier et all (2007, 
2008); with detailed references therein). Despite the reduction to 
a single wave form, the combined information of the GW ampli- 
tude and the location of the narrow peak of the GW spectral en- 
ergy density in frequency space contains information that makes 
it still possible to constrain progenitor- and postbounce rotation, 
but can barel y distinguish betwee n different finite-temperature 
EoS (see e.g. ID immelm eier et al] d2008l). who u sed the EoS of 
lLattimer & Swestvl dl99ll) and lShen et alJdl998bl) ). 

Gravitational waves fr om magneto- r otationa l collapse were 
considered in detail by IKotake et al] d2004bl) : IShibata et aD 



(2006j) and lObergaulin ger et all d2006h . As main differences 
compared to simulations which do not include magnetic fields 
it was found by these groups that only in the very special case 
of precollapse fields as s trong as > 10 12 G the overall dynam- 
ics can be influenced. As iHeger et alJ §005) have argued, such 
strong fields are unlikely to occur in standard core-collapse su- 
pernova progenitors. The GW amplitude is then affected by (i) 
time-dependent magnetic fields, which contribute considerably 
to the overall energy density, and (ii) by bipolar magnetohydro- 
dynamic (MHD) jet out flows which give rise to a so -called type 
IV signal with memory (Obergaulinger et al. 2006). Physically, 
such a memory effect in the GW signal arises from the temporal 
history of asymmetric matter outflow , leaving behind a constant 
offset in the amplitude dThornell989l) . 

Recently it has been argued through numerical simulations 
of equilibrium neutron star models or full core-collapse sim- 
ulations that PNSs with a high degree of differential rota- 
tion can be subject to nonaxiymmetric rotational instabilities at 
low p values (= T/\W\, ratio of rotational to gravit ational en- 
ergy) , leading to strong narrow-band GW emission dSaiio et alJ 
2003 : IWatts et alJl200a btt et al.ll200"l ISaiio & Yoshida) |2006: 
Ou & Tohlinell20q6t lOtt et al Jl2007bt ICerda-Duran et al.l l2007: 
btt et al.l l2007at IScheidegger et al.l 120081) . However little is 
known about the true nature of the instability at present. Previous 
work has failed to establish an analytical instability criterion, and 
the dependence of the instability on PNS rotation rate and de gree 
of differ ential rotation is still unclear, as it was pointed out bv lOttl 
d2009bh . 

Current state -of-the-art stellar evolution calculations 
(IHeger et al.ll2005l) tell us that iron cores of stars generally lose 
most of their angular momentum during their evolution due to 
magnetic torques. Therefore, they cannot become subject to 
strong rotationally-induced aspherities. However, anisotropic 
neutrino emission, convectio n and standin g accre tion shock 
instability (SASI, see e.g. iBlondin et all d2003l) ) - driven 
deviations from spherical symmetry which can lead to the 
emission of GWs of sizable amplitudes are likely to occur 
inside the PNS and the p ost-shock gain region and last for 
proba b ly hundreds o f ms dMiiller & Jankal 1 1997t iMiiller et al 
2004; Kotake etal J 120071: iMarek et alJ l2009n Murphv et al 
2009; Kotake et al.1 l2009allbir Although there is qualitative 
consensus among the different core-collapse supernova groups 
about the aforementioned features being emitters of stochastic 
broad-band signals, their detailed quantitative character still 
remains quite uncertain, since self-consistent 3D simulations 
with proper long-term neutrino transport were not carried out so 
far, but would in principle be required. The most recent, very 
elaborate 2D simu lations in that context were performed by 
IMarek et al.l ((2009). Their numerical setup includes long-term 
multi-flavor neutrino transport, an effective relativistic potential 
and two different EoS. The first 100ms of after bounce they 
observed GWs from early prompt postbounce convection, 
peaking somewhat around or below 100Hz, depending on the 
employed EoS. After this early episode, the GW emission in 
their models is dominated by a growing negative amplitude 
related to anisotropic neutrino emission at frequencies < 200Hz, 
while the GW signal associated with non-radial mass motions 
stems from density regimes 10 11 - 10 13 gcirT 3 and peaks in the 
frequency range of 300 - 800Hz, highly sensitive to the nuclear 
EoS. 

In the core-collapse scenario, PN S pulsations can provide 
another mechanism for GW emission. lOtt et al] (12006a) po inted 
out th at in the context of the acoustic mechanism (Burrows et al. 
2006) excited core g-mode oscillations might emit very strong 
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GWs. The GW signal is due to the nonlinear quadrupole compo- 
nents of t he pulsations tha t, at least initially, are of I = 1 g-mode 
character dOtt et al.ll2006al) . 

For recent reviews on GWs from cor e-co ll apse s uper- 
novae with comple te lists of references see lOttl (l2009bl) and 
iKotake et alJ (120061) . 

In this paper, we present the gravitational wave analysis of a 
comprehensive set of three-dimensional MHD core-collapse su- 
pernova simulations in order to investigate the dependencies of 
the resulting GW signal on the progenitor's initial conditions. 
Our calculations encompass presupernova models from stellar 
evolution calculations, a finite-temperature nuclear EoS and a 
computationally efficient treatment of the deleptonisation and 
neutrino emission during the collapse. General relativistic cor- 
rections to the spherically symmetric Newtonian gravitational 
potential are taken into account. Moreover, several models in- 
corporate long-term neutrino physics by means of a leakage 
scheme. As for the progenitor configuration, we systematically 
consider not only variations in the precollapse rotation rate, but 
also in the nuclear EoS and the magnetic field topology. This 
study extends the work of IScheidegger et all (120081) . who inves- 
tigated only two models with similar input physics. In this way, 
we carried out the so far largest parameter study of 3D MHD 
stellar collapse with respect to the prediction of GWs. 

The structure of the paper is as follows. In section 2 we 
briefly describe the numerical methods, input physics and pro- 
genitor configurations applied in our core-collapse simulations. 
Furthermore, we describe how we extract GWs from our model 
set. In sec. 3 we discuss results of 25 three-dimensional MHD 
simulations with respect to the GW signal from convection, rota- 
tional core bounce, nonaxisymmetric rotational instabilities and 
very strong magnetic fields. In sec. 4, we summarise and present 
conclusions. 



2. Numerical Methods 

2. 1 . The 3D MHD code and its input physics 

The algorithm used to solve the time-dependent, Newtonian 
MHD equations in our 3D simulations is based on a sim - 
ple and fast cosmologic al MHD code of iPen et al] J2003); 
Liebendo rfer et al.l d2006l) . which has been parallelised with a 
hybrid combination of MPI and openMP and improved and 
adapt ed to the requireme nts of core-collapse supernova simula- 
tions (Ka ppeli et al.ll2009t) . The 3D computational domain con- 
sists of a central cube of 600 3 km 3 volume, treated in equidis- 
tant Cartesian coordinate s with a grid spacing of 1km. It is, 
as explained in detail in Scheidegger et alJ d2008f) . embedded 
in a larger spherically symmetric computational domain that 
is treated by the time-i mplicit hydrodynamics code 'Agile' 
(Liebendorfer et al. 2002). With this set up, the code scales nicely 
to at least 8000 parallel processes (Kapp eli et alJ |200^). The 
ideal MHD equations read 



8E 

at 



dp 

■ (vpv - bb) + VP 

at 

+ V • \(E + P)\ - b • (v • b)] 

^ - V x (v x b) 
at 





-pVO 
-pv • 




(1) 

(2) 
(3) 
(4) 



expressing the conservation of mass, momentum, energy and 
magnetic flux, respectively, p stands for mass density, v is the 
velocity vector and E — pe + | v 2 + — the total energy (the sum 
of internal, kinetic and magnetic energy). The magnetic field is 
given by B — VAnb and the total pressure by P — p + y (the sum 
of gas and magnetic pressure). The MHD equations are evolved 
from a condition 



Vb = 0. 



(5) 



We e mploy the constrained transport method dEvans & Hawlevl 
1988) to guarantee the divergence-free time evolution of the 
magnetic field. The right hand side of Eq. [2] and Eq. [3] take 
into account the effect of gravitational forces on the magneto- 
hydrodynamical variables. The gravitational potential <t> obeys 
the Poisson equation 



A® = 4nGp 



(6) 



We only implement the monopole term of the gravitational po- 
tential by a spherically symmet ric mass integration that includes 
general relativistic corrections (Mare k et alj|2006i) . In this ap- 
proach, the Newtonian gravitational potential <£ is replaced by 
an effe ctive potential O c yy ('case A' potential, see iMarek et al.l 
(200f|) 



4> e //(r) 



r — 

Jr ~ 



nieff ,i 

1 + r' 3 -(p + p v ) 



An 



1 I p + e + p\ 



(7) 



with p being the gas pressure, p v the neutrino pressure and e the 
internal energy density. The effective mass is given by 

m eff (r) = An f dr V 2 (p + e + E) , (8) 
Jo 

where E is the neutrino energy. The metric function F is given 
by 



r= . i + v 2 - 



2m, 



eff 



(9) 



where v r is the radial fluid velocity. Note that lMiiller et all (2008) 
recently proposed an effective relativistic gravitational potential 
for rapidly rotating configurations. However, since it would have 
added another degree of freedom to our parameter set, we did not 
consider it in this study. 

The system of the MHD equations must be closed by a 
finite-temperature EoS. In core-collapse supernova simulations, 
an EoS has to handle several different regimes. For tempera- 
tures below 0.5MeV, the presence of nuclei and time depen- 
dent nuclear processes dominate the internal energy evolution. 
For simplicity, an ideal gas of Si-nuclei is assumed, which, for 
our GW study, is sufficiently accurate to describe the low bary- 
onic pressure contribution at low densities in the outer core. 
At higher temperatures, where matter is in nuclear statistical 
equilibrium, we employ two alternative EoS: The Lattimer & 
S westy EoS C attimer & Swestvl d!991l) . LS EoS) and the one 
bv lShen et alJ dl998al) . The first EoS is based on a phenomeno- 
logical compressible liquid drop model; it also includes surface 
effects as well as electron-positron and photon contributions. 
The LS EoS assumes a nuclear symmetry energy of 29.3MeV 
and we will perform simulations with three choices of the nu- 
cle ar compressibility modulus K (180,220,375MeV) as provided 
by lLattimer & Swestvl dl99ll) . Since variations in K affect the 
stiffness of the nuclear component of the EoS, this enables 
us to probe the effects of variations in stiffness while keep- 
ing the general EoS model fixed. The pure baryon EoS from 
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IShen et al.1 (Il998ah has a compressibility of K=281MeV and a 
symmetry energy of 36.9MeV. It is based on relativistic mean 
field theory and the Thomas-Fermi approximation. For matter 
in non-NSE (T< 0.44MeV), the Shen EoS is coupled to the 
baryon ic a nd electron-positron EoS g iven in iTimmes & Arnettl 
(1999j) and ITimmes & Swestvl d2000l) . It employs an ideal gas 
for the nuclei and additionally includes contributions from ion- 
ion-correlations and photons. 

2.2. The treatment of neutrino physics 

The treatment of neutrino physics is an essential ingredien t 
of core-collapse supernova simulations (Mezzacappa 2005). 
Multidimensional core-collapse supernova simulations therefore 
must rely on more or less severe approximations of the neutrino 
physics. In state-of-the-art 2D simulations, one way in which 
Boltzmann tra nsport is app roximated is the so-called 'ray -by-ray 
plus scheme' ( Bur as et al.l 12006!) . It solves the full transport in 
separate ID angular segments, where the neighbouring rays are 
coupled. Other groups rely o n multi-group flux limit ed diffusion 
(MGFLD, lOtt et al.l d2008b ; ISwestv & Mvral (l2009h ). MGFLD 
treats all neutrinos in seperate energy groups, drops the momen- 
tum space angular dependence of the radiation field and evolves 
the zeroth mom ent of the s pecific intensity instead of the specific 
intensity itself. lOtt et all (120081) also compared MGFLD with 
angle-dependent (i.e. partly Boltzmann) transport in 2D. In three 
dimensions, simulat ions have been perf ormed using 'grey' flux- 
limited diffusion dFrver & Warrenll2004l) . which oversimplyfies 
the important neutrino spectrum. It is important to resolve the 
neutrino spectrum, since the charge-current interaction rates go 
with the square of the neutrino energy. 

In our 3D MHD simulations, we app ly a parametrised delep- 
tonisation scheme dLiebendorferl 120051) . Detailed spherically- 
symmetric collapse calculations with Boltzmann neutrino trans- 
port show that the electron fraction Y e in different layers in the 
homologously collapsing core follow a similar deleptonisation 
trajectory with density. The local electron fraction of a fluid ele- 
ment can be parametrised as a function of density Y e (p). For our 
3D simulations, we apply tabulated Y e profiles which were ob- 
tained from detailed general relativistic, spherically symmetric 
three-flavour Boltzmann neutrino transport. Exemplary profiles 
are shown in Fig.Q] The offset between the Y e profiles from the 
LS- and the Shen EoS is consistent with the different asymme- 
try energies of the two EoS. This in turn is reflected in different 
neutrino reaction rates and thus a di fferent Y e at a given de nsity. 
These effects have b een discussed in Sumivoshi et al. (2008]) and 
iFischer et al.l d2009l) for massive progenitor stars in the range of 
4O-5OM . However, this parametrisation scheme is only valid 
until a few milliseconds after bounce, sin ce it cannot accoun t 
for the neutronisation burst, a s explained in Liebendorferl d2005l) 
and lScheidegger et all (120081) . 

In order to give an estimate of how long the above approxima- 
tion is able to predict quantitatively reliable postbounce (pb) GW 
signals for times t > 5ms after bounce, we carried out three rep- 
resentative simulations (R1E1CA L , R3E1AC L and R4E1FC L , 
see Tab. [TJ that incorporate a postbounce treatment for neutrino 
transport and deleptonisation. The scheme we apply is based 
on a partial implement ation of the isotropic diff usion source 
approximation (IDSA; iLiebendorfer et al.l d2009l) '). The IDSA 
splits the distribution function / of the neutrinos into two com- 
ponents, a trapped component/' and a streaming component/ 5 , 
representing neutrinos of a given species and energy which find 
the local zone opaque or transparent, respectively. The total dis- 
tribution function is the sum of the two components, / = f +f s . 




P [g/cm J ] 



Fig. 1. Electron fraction Y e as a function of density, obtained 
from detailed general relativistic, spherically symmetric three- 
flavour Boltzmann neutrino transport. Th e calculations were 
carried out using a 15M Q progenitor from Wo oslev & Weavei 
( 19951) an d a finite temperatur e EoS of either lLattimer & Swestv 
(19911) or lShen et all dl998aj). In the figure, we denote the LS 
EoS version with a compressibility of 180MeV as LSI 80, the 
one carried out with 375MeV as LS375. The Shen EoS is named 
as such. 



The two components are evolved using separate numerical tech- 
niques, coupled by a diffusion source term. The trapped com- 
ponent transports neutrinos by diffusion to adjacent fluid ele- 
ments, while in this paper the streaming component is discarded 
(/■' = 0), implying that these neutrinos are lost from the simula- 
tion immediately. The part of the IDSA for trapped neutrinos is 
implemented in three dimensions, including both electron neu- 
trinos and electron anti-neutrinos (Whitehouse & Liebendorfer 
2010, in preparation). The use of this partial implementation of 
the IDSA enables us to capture the neutrino burst and to continue 
our simulations into the postbounce regime, as shown in Fig. [2] 
The leakage scheme is switched on at core bounce, when the 
central density of the core reaches its first maximum value. At 
present, we neglect the emission of /j and r neutrinos and their 
antineutrinos, which in principle play an important role in the 
cooling of the PNS. However, since our leakage schemes ne- 
glects any absorption of transported neutrinos, it already overes- 
timates the cooling of the PNS, even without the treatment of the 
/j. and t neutrinos. This is shown in [3] 

2.3. Initial model configurations 

We construct the initial conditions of our simulations by a para- 
metric approach. All our models are launched from a 15M Q 
progenitor star from the stellar evolutionary calculations of 
IWooslev & Weaverl dl995l) . Angular momentum was added to 
the presupernova model according to a shell-type rotation law 
(lEriguchi & Muller|[T98l 



Q(r) = Ru 



A 2 + r 2 



(10) 



where we define r = ^x 2 + y 2 + z 2 as spherical radius, and 
(x,y,z) as Cartesian coordinates. The constant A is the degree 
of differential rotation that controls the steepness of the angu- 
lar velocity profile, Q,- e the initial central rotation rate, and r is 
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0.4 0.6 

Mass [Ms] 



Fig. 2. Comparison of the Y e profiles of the almost non-rotating 
3D models R1E1CA and R1E1CA L with the spherically sym- 
metric model G15, which is based on general relativistic three- 
flavour neutrino Boltzmann transport ( Lieb endorfer et ai] |2005) 
as function of the enclosed mass at 5ms before and after bounce. 
Note that while the simulation R1E1CA, which includes only the 
neutrino parametrisation scheme for the collapse phase, cannot 
model the neutrino burst, the 'leakage' -IDS A model RIEICAl 
can. 



the distance from the origin. In our whole model set we choose 
A = 500km as degree of differential rotation. 

In order to guarantee a divergence-free initial state, the initial 
magnetic field configuration was set up employing its definition 
via the vector potential A. The components of A we chose are 



A = (-^y,^x,J^B tor ) 



(11) 



In order to mimic a dipole-like field, we scale the vector potential 
with density according to 



Pref 



(12) 



Finally, the magnetic field is derived from this vector potential 

B = VxA. (13) 

The initial toroidal- and poloidal components of the magnetic 
field are spe cified at a reference density of p re f = 5 x 10 7 gcm~ 3 
according to lHeger et al.l d2005l) . 

We compute a total of 25 models, changing the combination 
of total angular momentum, the EoS, and toroidal- and poloidal 
magnetic fields. The model parameters are summarised in table 
Q] The models are named after the combinations of initial central 
rotation rate, the EoS, and toroidal- and poloidal magnetic fields. 
The first two letters of the model name represent the initial cen- 
tral rotation rate Q c>! [rads -1 ] according to 






0.3 


71 


2tt 


3tt 


An 


R0 


Rl 


R2 


R3 


R4 


R5 



the second two letters stand for the applied EoS 



LS 


LS 


LS 


Shen 


(K = 180MeV) 


(K = 220MeV) 


(K =325 MeV) 




El 


E2 


E3 


ST 



while the last two letters assign the order of magnitude of the 



toroidal-. 

10 6 



and poloidal field strength [G] according to 



10 7 



10" 



10" 



10' 



D 

The final subscript L signs simulations which were carried out 
with the leakage scheme. Some of the values which we adopt as 
rotation r ate and magnetic fi elds correspond to the values sug- 
gested in Hegeret al. (2005). However, we point out that some 
of the initial rotation rates (Q ty > 7rrads _I ) and magnetic fields 
(B > 10 10 G) are larger and stronger compared to current pre- 
dictions from stellar evolution calculations. However, we com- 
puted these models in order to cover a wide parameter space for 
our three-dimensional models with neutrino transport approxi- 
mations. 



2.4. Gravitational Wave extraction 

The two independent polarisations of the dimensionless gravita- 
tional wave field h^ v = g MV - rj^ v in the transverse traceless gauge 
are given by 

hJj r (X,t) = ^(A + e + +A x e x ). (14) 

The spatial indices i, j run from 1 to 3 (Mis ner et alJll973"h . R is 
the distance from the source to the observer and the unit polari- 
sation tensors e+ and e x in spherical coordinates are represented 

as 

e+ = e <g> e g - e$ ® e$ (15) 
e x = e e ®e 4 , + e 4 ,®e e . (16) 



In th e slow-motion limit dMisner et al.l 1 1 9731: iFinn & Evans! 
1 1990b the amplitudes A + and A x are linear combinations of 
the second time derivative of the transverse traceless mass 
quadrupole tensor tTT . In Cartesian coordinates, the quadrupole 



tensor is expressed as 



-6ijix\ 



A) 



(17) 



For convenience we evaluate below the GW amplitudes along 
the polar axis (9 — <p = 0, denoted as subscript I) 



■JTT 



■JTT 

yy 



A + i .., 

and in the equatorial plane (0 - |, (f> = 0, denoted as II) 



A.n = fj - t 



■JTT 



ixll 



-2f T 



(18) 
(19) 

(20) 
(21) 



Since a direct e valuation of Eq . (fT71) i s numerically problematic, 
as discussed in IFinn & Evansl (Il990l) . we apply alternative re- 
formulations of the standard quadrupole formula, in which one 
or both time derivatives are replaced by hydrodynamic variables 
using the continuity and momentum eq uations. In the first m o- 
ment of momentum density formula (see lFinn & Evansl JT990)). 
the first time derivative of the quadrupole moment yields 



; TT 



4/ 



dVp 



VjXj + vjXi - -6ij (viJCi + v 2 x 2 + v 3 x 3 ) 



.(22) 



In the stress formulation iBlanchet et al.1 (1 19901) . f?7 is given by 
r iV = ^ / dVp(2 Vi vj - Xi dj® eff - Xjd i O eff ) TT , (23) 
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Table 1. Summary of initial conditions 



Model 


ii Cj i |rads J 


Pi 




Pb 




EoS 


.-. r S 


-] 




B po ij[G] 


E m /\W\, 




E m /\W\ f 


R0E1CA 















El 


4.39 x 


nl4 
10 


5 x 10 9 


1 x 10 6 


1.8 x 10" 


-9 


3.8 x 10 


R0E3CA 















E3 


4.17 x 


r\14 

10 


5 x 10 9 


1 x 10 6 


1.8 x 10" 


-9 


3.1 x 10~ 8 


ROSTCA 















ST 


3.38 x 


r>l4 

L0 14 


5 x 10 9 


1 x 10 6 


1.7 x 10" 


-9 


3.3 x 10~ 8 


R1E1CA 


0.3 


0.59 x 10 


_5 


1.7 x 10 


-4 


El 


4.53 x 


10 14 


5 x 10 9 


1 x 10 6 


1.7 x 10" 


-9 


3.4 x 10 


R1E3CA 


0.3 


0.59 x 10" 


_5 


1.7 x 10 


-4 


E3 


4.17 x 


10 14 


5 x 10 


1 x 10 6 


1.7 x 10" 


-9 


3.6 x 10~ 8 


R1E1DB 


0.3 


0.59 x 10 


-5 


1.7 x 10 


-4 


El 


4.53 x 


nl4 

10 


5 x 10 10 


1 x 10 7 


1.7 x 10' 


-7 


3.0 x 10 


RISTCA 


0.3 


0.59 x 10 


-5 


1.7 x 10" 


-4 


ST 


3.38 x 


L0 


5 x 10 9 


1 x 10 6 


1.7 x 10" 


-9 


3.2 x 10~ 8 


R1E1CA L 


0.3 


0.59 x 10 


_5 


1.8 x 10" 


_4 


El 


4.38 x 


to 14 


5 x 10 9 


1 x 10 6 


6.1 x 10" 


-9 


1.2 x 10~ 7 


R2E1AC 


n 


f\ s~ A , , 1 f\- 

0.64 x 10 


_3 


1.6 x 10 


-2 


El 


4.27 x 


10 14 


1 x 10 6 


1 x 10 9 


6.5 x 10" 


-9 


1.27 x 10~ 5 


R2E3AC 


n 


0.64 x 10 


_3 


1.6 x 10" 


-2 


E3 


4.00 x 


nl4 

10 14 


1 x 10 s 


5 x 10 9 


6.4 x 10" 


-9 


5.1 x 10~ s 


R2STAC 


7X 


0.64 x 10' 


_3 


1.6 x 10 


-2 


ST 


3.30 x 


n 14 

10' 4 


1 x 10 6 


5 x 10 J 


6.4 x 10" 


-9 


3.6 x 10~ 6 


R3E1AC 


2n 


0.26 x 10' 




5.2 x 10" 


-2 


El 


3.80 x 


10 14 


1 x 10 6 


5 x 10 9 


8.0 x 10" 


-9 


6.7 x 10~ 6 


R3E2AC 


2n 


0.26 x 10' 




5.1 x 10 


-2 


E2 


3.65 x 


10 14 


1 x 10 6 


5 x 10 9 


8.9 x 10" 


-9 


7.2 x 10~ 6 


R3E3AC 


2n 


0.26 x 10' 


_2 


5.1 x 10 


-2 


E3 


3.64 x 


10 14 


1 x 10 6 


5 x 10 9 


9.0 x 10" 


-9 


5.8 x 10 


R3STAC 


2n 


0.26 x 10' 


_2 


5.1 x 10" 


2 


ST 


3.01 x 


nl4 
10 14 


1 x 10 6 


5 x 10 9 


8.9 x 10" 


-9 


4.7 x 10 


R3E1CA 


2n 


0.26 x 10' 


_2 


5.2 x 10" 


-2 


El 


3.82 x 


10 14 


5 x I0 J 


1 x 10 6 


3.1 x 10 


-9 


1.4 x 10" 8 


R3E1DB 


2n 


0.26 x 10' 


_T 


5.2 x 10" 


-2 


El 


3.81 x 


n]4 

10 


5 x 10 10 


1 x 10 7 


3.1 x 10" 


-7 


1.6 x 10~ 6 


R3E1AC L 


2n 


0.26 x 10' 


_2 


5.1 x 10" 


-2 


El 


3.65 x 


nl4 

10 14 


1 x 10 6 


5 x 10 9 


8.7 x 10" 


-9 


8.6 x 10 -5 


R4E1AC 


3n 


0.57 x 10' 


_2 


8.6 x 10" 


_2 


El 


3.22 x 


10 14 


1 x 10 6 


5 x \W 


1.2 x 10" 


-8 


1.3 x 10~ 5 


R4STAC 


3n 


0.57 x 10' 


_2 


9.0 x 10" 


2 


ST 


2.69 x 


Q14 


1 x 10 6 


5 x 10 J 


1.5 x 10" 


-8 


1.4 x 10~ 5 


R4E1EC 


3n 


0.57 x 10" 


_2 


8.6 x 10" 


-2 


El 


3.22 x 


nl4 

10 


1 x 10" 


5 x 10 9 


2.4 x 10' 


-8 


1.3 x 10" 5 


R4E1FC 


3n 


0.57 x 10' 


_2 


8.6 x 10 


-2 


El 


3.19 x 


to 14 


1 x 10 12 


5 x 10 9 


1.3 x 10' 


-4 


1.3 x 10~ 4 


R4E1FQ 


3n 


0.57 x 10' 




8.7 x 10 


-2 


El 


3.14 x 


10 14 


1 x 10 12 


5 x 10 9 


1.3 x 10' 


-4 


2.0 x 10- 4 


R4E1CF 


3n 


0.57 x 10' 


_T 


8.2 x 10 


-2 


El 


3.22 x 


LO 14 


5 x 10 9 


1 x 10 12 


5.8 x 10 


-4 


5.9 x 10" 3 


R5E1AC 


An 


1.02 x 10- 




10.2 x 10 


_2 


El 


2.47 x 


10 14 


1 x 10 6 


5 x 10 9 


1.3 X 10' 


-8 


1.2 x 10- 5 



" The subscript i stands for initial, b for bounce, while / stands for final, pi, is the maximum central density at the time of core bounce. B po ij 
and B torj abbreviate the initially imposed toroidal and poloidal magnetic fields, whereas E m /\W\ stands for the ratio of magnetic to gravitational 
energy. The subscript L denotes models which were carried out with the leakage scheme. 



where Q> e ff is the gravitational potential. We apply both ways of 
computin g the GW signal and co mpare them by computing the 
overlap (iGiacomazz o et al. 2009) of the resulting waveforms via 



0(h u h 2 ) 



(hx\h 2 ) 



(24) 



{h x \h 2 ) = 4% df 



f 

Jo 



C M/)W) 
s h (f) 



(25) 



where 5/,(v) [Hz -1 ] is the power spectral density of the strain 
noise from a given detector, e.g. LIGO, which were kindly pro- 
vided by Shoemaker (2007, private communication) The results 
show good agreement with numerical deviations of a few per- 
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Table 2. Overlap 



Model 


Q. c i [rads '] 


Overlap 


R1E1AC L 


0.3 


96% 


R3E1AC L 


2k 


92% 


R4E1AC L 


3n 


97% 



cent, as displayed in Table|2] There, we computed the overlap of 
GW trains from the representative models R1E1CA/,, R3E1CA£ 
and R4E1FCl, extracted by using Eq.l22land Eq.[23] We assume 
optimal orientation of detector and source as well as the GW to 
be emitted along the polar axis. Below, we generally apply Eq. 
[22] if not stated otherwise. Given our spherically-symmetric ef- 
fective GR approach to the solution of the Poisson equation, this 
expression is physically best motivated, since it does not depend 
on spatial derivatives of the gravitational potential. 

We also point out that the Newtonian quadrupole formalism 
to extract the gravitational radiation is not gauge invariant and 
only valid in the Newtonian slow-motio n limit. However, it was 
shown by Shi bata & Sekiguchl ([2003) that the above method 
seems to be sufficiently accurate compared to more elaborate 
techniques, as it preserves phase while being off in amplitude 
by ~ 10% in neutron star pulsations. The energy carried away 
by gravitational radiation can be calculated by the following ex- 
pression: 



Ecw = hf[jt{ IiJ ~^ ijI '')} dt 

15G J L ■« yy & 

l.xx^yy ~ Ixx^zz ~ ^yy^zz 

+ 3(4 + 4 + ^)] ' 

where 7, ; = ?.J r . The equivalent frequency integral yields 



(26) 



4c 3 
15G 



2 + p +p 

yy zz 



'xx*yy Ixx^zz lyy^zz 



3 (p +/ 2 + 

\ xy xz 



(27) 



where /,y( v ) * s tr, e Fourier transform of the quadrupole amplitude 

In order to calculate the contribution to the GW signal due to 
magnetic stresses, we generalised Eq. d23l , taking into account 
contributi ons from the magn etic field. Following the deriva- 
tions of iKotake elatl d2004al) and lObergaulinger et ail (§006), 
the Cartesian quadrupole gravitational wave amplitude in the 
MHD case yields 

W = ^f dV [ 2 fu ~ P (ti^eff + XjdPeff)]"' , (28) 

where fy = pv(Vj—bib j, and = B-J yf^n. Below, the GW source 
is assumed to be located at the Galactic centre at R = lOkpc. 
We estimate the signa l-to-noise ratios (SNR) for o ptimal filtering 
searches according to Flanagan & Hughes ( 1998): 



iMv)i 2 + iMv)r 

S h (v) 



dv , 



(29) 



where h(v) is the Fourier transform of the gravitational wave am- 
plitude, and Si,(v) is the power spectral density of strain noise 
in the detector. We assume optimal orientation of detector and 
source, and the Fourier spectra were normalised according to 
Parseval's theorem. 



3. Results 

In the subsequent four subsections we will discuss the GW sig- 
nature of our 25 models. While presenting the resulting GW pat- 
terns, we will pay special attention to possible imprints of dif- 
ferent finite temperature EoS, rotation rates, nonaxisymmetric 
instabilities, magnetic fields and a postbounce leakage scheme 
on the predicted 3D GW signals. The models' initial conditions 
and similar relevant quantities are summarised in Table[TJ whilst 
the GW data is listed in Table [3] 



3.1. Non- or slowly rotating core collapse 
General remarks 

Non- and slowly rotating progenitors (Q tv = . . .0.3rads~' in 
our model set) all undergo quasi-spherically symmetric core col- 
lapse. As the emission of GWs intrinsically depends on dynami- 
cal processes that deviate from spherical symmetry, the collapse 
phase in our models (that neglect inhomogenities in the progen- 
itor star) therefore does not provide any kind of signal, as shown 
in Fig.[3]for t-tt, < 0. However, subsequent pressure-dominated 
core bounce, where the collapse is halted due to the stiffening 
of the EoS at nuclear density p nuc m 2 x 10 14 gcirT 3 , launches a 
shock wave that plows through the infalling material, leaving be- 
hind a negative entropy gradient that induces so-called 'prompt ' 
connective activity ( e .g.lMiiller et al.l(l2004t):lBuras et al](| 2006): 
iDimmelmeier et al.l (l2008l) : iMarek et al l (I2009I) : lOttl (l2009b)). 
The GW burst which is accompanied by such aspherities starts 
several ms after bounce when convective overturn starts to be 
effective. The criterion for convective instability (the 'Ledoux 
condition') is generally expressed as (lLandau & Lif shit3 1 1 9 5 9b 
IWilson & Mavle 198§ 



dp_ 

dY e 



<" IP,s 



dYe 

dr 



P,Y. 



£)><>• 



(30) 



where Y e , p, r and s are electron fraction, density, stellar 
radius and entropy per baryon respectively. Thermodynamic 
consistency requires (dp/ds)y e ,p < 0, which implies that a 
negative entropy gradient always acts in a destabilising man- 
ner. Additionally, the neutronisation burst, occurring some ~ 
5ms after bounce, causes a negative lepton gradient at the 
ed ge of the PNS w h ich fu rt her drives c o nvecti o n (cf. Fig. 13 
of^S umivoshi etal.1 (120041) . iBuras et all (120061) : iDessart et all 
(120061) 1. Note that convection is be weakened in the rotational 
plane by positive specific angu lar momentum gradients in rotat- 
ing cores (Endal & Sofial ll978l) . 



Models without deleptonisation in the postbounce phase: 
Effects of the EoS and magnetic fields on the GW signature 

During the early postbounce stage (t — % < 20 - 30ms), prompt 
convective motion is predominantly driven by a negat i ve en- 
tropy gra dient, as pointed out e.g. by IScheidegg er et al.1 (12 009) 
see als o rDimmel meier et al.l (l2008h : IMarek et al.l (120091) : lOttl 
2009b)). The GW burst to be associated with prompt convec- 
tion sets in ~ 6ms after bounce in models based on the LS EoS, 
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Table 3. Summary of GW related quantities 



Model 


f/[ms] 




dir. 






l-^+./^.mfl.vl 


\Ax,pk,max\ 


MHz] 


Mm 


R0E1CA 


130.8 


2.15 x 10-" 


I 
II 


3 
2 


< 1 

< 1 


3 
5 


< 1 
1 


- 


- 


R0E3CA 


103.9 


5.72 x 10-" 


I 
II 


2 
2 


< 1 

< 1 


10 

8 


< 1 

< 1 


- 


- 


ROSTCA 


70.3 


1.35 x 10-" 


I 

n 


3 
1 


< 1 

< 1 


4 

3 


< 1 

< 1 


- 


- 


R1E1CA 


112.8 


1.36 x lO" 10 


i 
ii 


3 
3 


< 1 
<1 


12 
15 


6 
1 


- 


- 


R1E3CA 


130.4 


1.43 x 10" 10 


i 
ii 


2 
3 


< 1 

< 1 


10 
17 


5 

< 1 


- 


- 


R1E1DB 


112.8 


1.24 x 10-'° 


i 
ii 


3 
3 


< 1 

< 1 


12 
15 


6 
2 


- 


- 


R1STCA 


45.8 


2.01 x 10-" 


i 
ii 


2 
2 


<1 
<1 


3 
6 


1 

<1 


- 


- 


RlElCAz, 


92.9 


1.04 x lO" 10 


i 
ii 


3 
2 


< 1 
<1 


4 
10 


2 
1 


- 


- 


R2E1AC 


127.3 


5.52 x 10- 9 


i 
ii 


2 
105 


<1 


6 
1 


5 

3 


841 


- 


R2E3AC 


106.4 


5.31 x lO" 5 


i 

n 


2 
104 


<1 


4 
1 


4 
< 1 


803 


- 


R2STAC 


64.0 


7.62 x 10" 9 


i 

ii 


3 

133 


<1 


6 
2 


6 
1 


680 


- 


R3E1AC 


62.5 


5.58 x 10- 8 


i 
ii 


2 

393 


< 1 


15 
10 


16 

2 


880 


725 


R3E2AC 


45.5 


6.53 x 10~ 8 


i 

n 


2 

409 


< 1 


11 

5 


13 

3 


882 


- 


R3E3AC 


57.5 


6.44 x 10~ 8 


i 

ii 


2 

409 


< 1 


6 
6 


9 

3 


891 


- 


R3STAC 


50.2 


1.05 x 10~ 7 


i 
ii 


5 

526 


< 1 


9 
5 


11 

3 


854 


935 


R3E1CA 


69.4 


8.05 x 10~ 8 


i 

n 


3 

426 


< 1 


9 
5 


9 
3 


897 


- 


R3E1DB 


62.7 


7.76 x 10~ 8 


i 

ii 


8 

425 


< 1 


13 
8 


11 

3 


886 


- 


R3E1AC L 


196.7 


2.14 x 10- 7 


i 
ii 


1 

437 


< 1 


136 
70 


136 
9 


909 


909 


R4E1AC 


98.7 


7.74 x 10~ 8 


i 

n 


1 

512 


< 1 

< 1 


37 
20 


37 
6 


385 


662 


R4STAC 


67.2 


1.91 x 10~ 7 


i 
ii 


2 

536 


< 1 


102 
53 


109 

3 


396 


902 


R4E1EC 


100.8 


7.51 x 10~ 8 


i 
ii 


1 

492 


< 1 


25 
15 


20 
7 


866 


611 


R4E1FC 


80.1 


7.29 x 10~ 8 


i 

n 


< 1 
516 


< 1 

< 1 


16 
9 


19 

3 


859 


828 


R4E1FC L 


97.6 


3.42 x 10~ 7 


i 
ii 


2 

536 


< 1 


316 
157 


298 
6 


485 


673 


R4E1CF 


19.1 


6.50 x 10- 8 


i 
ii 


1 

518 


< 1 


4 

ring-down 


4 

2 


370 




R5E1AC 


93.2 


1.20 x 10- 8 


i 

n 


1 

238 


< 1 


18 
10 


18 
11 


317 


727 



" tf is the time after core bounce when the simulation was stopped. Egw is the total energy released in gravitational radiation. We present 
the maximum amplitudes at different stages of their time-evolution in polar (I) and equatorial (II) direction. The subscripts b and pb stand for 
bounce and postbounce. fi, denotes the peak frequency of the GW burst at bounce, while f TW stands for the spectral peak from the narrow band 
emission caused by a low T/\W\ instability. The electronic wave forms of all models are available upon request. 



generally ~ 1 — 3ms before the same feature occurs in the cor- 
responding simulations using the Shen EoS, as indicated in Fig. 
[3] The reason for this behaviour is that the shock wave in the 
'Shen' -models carries more energy compared to those in mod- 
els using a LS EoS. Therefore, the shock wave stalls at slightly 
later times at larger radii (see Fig.|4]i, and the conditions for con- 
vective activity are delayed compared to the LS runs. Note that 



the convective overturn causes a smoothing of the negative en- 
tropy gradient. As a result, the GW amplitude quickly decays 
(f > 30ms after bounce) and is not revived during the later evo- 
lution of the models without deleptonisation in the postbounce 
phase, as displayed in Figs.[3]and[5] Simulations that incorporate 
the Shen EoS return up to a factor of 2 smaller maximum am- 
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Fig. 3. Time evolution of the GW amplitude A+i from the slowly 
rotating models R1STCA (full line) and R1E1CA (dashed line). 
Since convection is a stochastic process, the GW amplitudes are 
rather insensitive to the location of the observer. Hence, we dis- 
play only one representative polarisation. 




0.2 0.4 0.6 0.8 1 1.2 

Mass [Ms] 



Fig. 4. Spherically averaged density profiles from the slowly ro- 
tating models R1E1CA (red full line), R1E3CA (black dotted 
line) and R1STCA (black dashed line) at bounce. The second 
entropy time slice is chosen to be approximately at the onset of 
the GW signal from prompt convection. 

plitudes compared to their counterparts, as can be deduced from 
Fig. [3] and Tabled 

We also find that the GWs from simulations which were car- 
ried out with the stiff E3 show no significant deviations from 
those computed with El (see Tableland Fig . O . The minor de- 
viations of the GWs are entirely due to the stochastic nature of 
convection. The waveform spectra from the LS models cover a 
broad frequency band ranging from ~ 150 - 500Hz, as displayed 
in Fig. [6] The spectral peak of the Shen models is shifted some- 
what to lower frequencies, covering a range from ~ 150- 350Hz 
(see Fig. |6). When comparing the energy Eqw which is emitted 
during the first 30ms after bounce from the LS models to that of 
the corresponding Shen models, we find the latter models emit 
less energy (see Table |3). This discrepancy is due to the lower 
emission at higher frequency v > 350Hz in the Shen models and 
the fact that the emitted energy dEcw/dv is proportional to v 2 . 



Fig. 5. Time evolution of the GW amplitude A + n from the non- 
rotating models R0E1CA (dashed) and R0E3CA (full line). 

Moreover we find that already slow rotation (rotation rate Rl) 
leads to a deformation of the PNS. This can be quantified by 
considering e.g. a density cut at p — 10 n gcm -3 , which is just 
inside the convectively unstable region. For the slowly rotating 
model R1E1CA, this point is located at a radial distance of 80km 
along the polar axis at 20ms after bounce. However, in the equa- 
torial plane, the same position is reached one radial grid zone 
further out relative to the origin due to the action of centrifugal 
forces. The short time variation in the quadrupole due to rotation 
combines with that of prompt convection and together they lead 
to somewhat stronger GW emission in the slowly rotating case 
compared to the non-rotating model set. This effect is strongest 
for the GW amplitude A + n, which is the 'axisymmetric' (1=2, 
m=0) component of the wave field. Despite this feature, the fre- 
quency content of models which only differ in rotation rate, stays 
practically unaltered. 

We found the key controlling factors that govern the GW 
emission from prompt convection to be the i) radial location of 
the convectively unstable zones and ii) the related characteristic 
dynamical timescales involved, for which we use as rough es- 
timate tdyn ~ A r /c7E] (Ott 2009, private communication). Note 
that both i) and ii) are implicitly determined by the applied EoS, 
being responsible e.g. for the local speed of sound and the ra- 
dial PNS density profile. Our core-collapse simulations that use 
a version of the LS EoS (El or E3) show at core bounce maxi- 
mum central densities up to ~ 25% higher than the correspond- 
ing models that apply the Shen EoS (see TabQ}. Moreover, the 
LS-models possess a PNS which is more strongly condensed in 
central regions and has a steeper density gradient further out. 
The densities of e.g. R1E1CA and R1STCA intersect at ~ 0.5M Q 
(~ 8km), as displayed in Fig. [7] In addition, models using the El 
show somewhat higher central densities compared to the ones 
applying the stiffer E3 variant, and densities of e.g. R0E1CA 
and R0E3CA cross at ~ 7km. Nevertheless these two variants of 
the LS EoS differ only in their compressibility. The radial struc- 
tures resemble each other strongly, unlike to the models carried 
out using the Shen EoS, as shown in Fig. [7] Hence, the partic- 
ular similarity in the GW characteristics of LS runs (see Fig. [5]) 
is not unexpected since the region which is convectively unsta- 
ble is restricted to roughly the same radial- and density regimes 

1 ~ = 1/A r j c s (r)dr is the radially averaged sound speed of a con- 
vectively unstable layer with a radial extension of A,-. 
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- R0E1CA; observer position: (6=7i/2,(^=0)-axis 

- ROSTCA; observer position: (9=7i/2,(^=0)-axis 
-LIGO 



AdvLIGO 
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Fig. 6. Spectral energy distribution from the models R0E1CA 
(dashed line) and ROSTCA (full line) for a spectator in the equa- 
torial plane at a distance of lOkpc compared with the LIGO 
strain sensitivity (Shoemaker 2007, private communication) and 
the planned performance of Advanced LIGO. Optimal orienta- 
tion between source and detector is assumed. 
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Fig. 7. Spherically averaged density profiles from models 
R1E1CA (dotted line), R1E3CA (full line) and R1STCA 
(dashed line) at core bounce (f - t/, = 0). 



(R * 30-70km,p a 10 11 - 1.1 x 10 12 gem" 3 ) and therefore bound 
to the same dynamical timescale td yn and its corresponding fre- 
quency band. For the Shen models on the other side the region 
of negative entropy gradient, whose boundaries are radially con- 
strained by R x 60-90km, contains considerably less matter than 
its LS analogues. This explains the smaller GW amplitudes. A 
narrower density spread (p « 7 ■ 10 10 - 1.2 x 10 11 gem -3 ) in the 
unstable region leads to a more restricted peak frequency band 
at lower values (see Fig. |4j». 

Computing the SNRs we find that all the simulations dis- 
cussed above lie just below the detector limits of LIGO if we 
assume them to be located at a Galactic distance of lOkpc. Their 
single-detector optimal-orientation SNR is just a little above 
unity. However, for a successful detection, at least a SNR of 



7 to 8 is necessary. Note that the SNRs of models with differ- 
ent EoS do not differ much at this stage since all share a sim- 
ilar spectral energy distribution within the window of LIGO's 
maximum sensitivity. The current detector sensitivity does not 
allow for the detection of the high frequency tail of the LS mod- 
els. Note, however, that for planned future detectors such as the 
Advanced LIGO facility, things change dramatically. As a di- 
rect consequence, these new detectors would permit the distinc- 
tion between the prompt convection GW signal from the LS and 
Shen EoS, since the full spectral information would be available. 
However, we find it impossible to discriminate between the dif- 
ferent LS EoS variants. Hence, our simulations indicate that the 
GW signature depends more strongly on the asymmetry energy 
than the compressibility parameter of the EoS. 

These results are partly di f ferent from those previousl y pub- 
lished. Recently. iMarek et al.1 d2009l) as well as[Ott (2009g) re- 
ported to observe GW from prompt convection in state-of-the- 
art 2D simulations which were launched from similar initial 
conditions as ours, namely the same progenitor star and the 
soft variant of the LS EoS (K=180MeV). Whilst t he extracted 
GW a mplitudes from early prompt convection of IMarek et al.l 
(2009) (cf. their model M15LS-2D) are in rough agreement with 
our results, the spectrum of their wave train peaks at consid- 
erably lower frequencies, namely around about ~ 100Hz. We 
suppose that this discrepancy is due mainly to different radial 
locations of the unstable regions and t he consequen tly encom- 
passed amount of overturning matter. lOtj (1200 9b) computed 
two models with different resolution. While one of their mod- 
els (sl5WW95) in particular fits our results well in all char- 
acteristic GW features, namely the size of amplitudes, band of 
emission and the amount of emitted energy, the better resolved 
model (sl5WW95HR) showed that convection is much weaker 
due to less seed perturbations and hence the GW signal and 
the total amount of emitted energy considerably lower. We re- 
cently also tes ted this issue with a be tter resolved model (cf. 
model RIhr of Scheidegger et al. (2009), with a grid spacing of 
0.6km). This better resolved model showed considerably smaller 
seed perturbations around t — % ~ 0, as grid alignment effects 
are better suppressed at core bounce; hence prompt convection 
then is much weaker and a smaller GW amplitude (~ 50%) is 
emitted. However, better numerical resolution also leads to less 
numerical dissipation in th e system, which eases th e dynami- 
cal effects that follow. Thus. IScheideggeret al. (2009) found for 
~ 10 < t < 20ms considerably stronger GW emission from early 
prompt convection compared to the 1km resolved models. 

Finally, we find that imposing initial magnetic fie l ds ten 
times as strong as the values suggested in iHeger et al.l ([2005 ) 
does not influence the model dynamics and therefore the GW 
signal at all, as can be deduced from Table[3]by a comparison of 
models R1E1CA and R1E1DB. 



Model with deleptonisation in the postbounce phase 

Due to the absence of accurate postbounce neutrino transport, 
it is unclear how reliably the models discussed in the previous 
subsection predict the GW signals from the early postbounce 
period. In order to investigate this question, we carried out one 
computationally expensive simulation, model R1E1CA/,, that in- 
cludes the emission of neutrinos after bounce but neglects the 
neutrino heating, which becomes relevant at t - % > 50ms. The 
absence of the latter makes it impossible for this model to track 
the long-term postbounce neutrino-driven convection. However, 
in comparison with purely MHD models, the inclusion of neu- 
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Fig. 9. Model RlElCA^'s spherically averaged specific entropy 
(full line) and 10 x electron fraction Y e (dashed line) profiles are 
plotted versus the enclosed mass ~ 5ms after bounce. The radial 
position of simultaneous negativ entropy- and lepton gradients 
are marked by the vertical lines. 
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Fig. 10. Model RlElCA/,'s spectral energy distribution of the 
GW signals at a distance of lOkpc. Note that the spectrum inter- 
val from ~ 150 - 500Hz is caused by prompt convective activity, 
while the higher modes result from lepton gradient driven post- 
bounce PNS convection. 



trino leakage makes is possible to distinguish effects on the GW 
signal due to entropy and lepton-gradient driven convection. 

A detailed comparison of model RIEICA^ with its purely 
magnetohydrodynamical counterpart R1E1CA shows that both 
follow a similar dynamical behaviour until about 20ms after 
bounce. Asphericities leading to GW emission are predomi- 
nantly driven by entropy- and not lepton-induced convection 
in this supernova stage. Consequently, the wave trains emitted 
within this early period fit each other qualitatively (cf. Fig. [3] 
and Fig . [HJ . However, we find some quantitative deviations: The 
GWs of RIEICA/, reach lower maximum values (see Table [3]), 
as the 'Ledoux' unstable region encompasses less mass (Fig. 
[5}, and the presence of neutrino cooling leads to a more rapid 



smoothing of the entropy gradient compared to the models dis- 
cussed in the previous subsection. However, since the over- 
turning matter in the top layers of the PNS has the same ra- 
dial position, densities and dynamical timescales as in model 
R1E1CA, these models have similar GW spectra, peaking be- 
tween ~ 150 - 500Hz (cf. Fig. and Fig. [10j>. Therefore the 
physically simpler models still provide reasonably accurate GW 
predictions in frequency space until about 20ms after bounce, 
although the amplitudes are overestimated a few xl0% (cf. 
Figs. [3] and |8]l. Model RlElCA^'s later postbounce evolution 
(f > 20ms) differs strongly compared to its purely hydrody- 
namical counterpart R1E1CA. A negative radial lepton gradi- 
ent, caused by the neutronisation burst and subsequent delep- 
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Fig. 12. Summary of all model's dimensionless peak gravita- 
tional wave amplitude \h\ nwx j, at a distance of lOkpc and at core 
bounce versus the rotation rate /3h- While \h\ max ^ scales roughly 
linearly with /3b for rotation rates /3b < 10%, the growing cen- 
trifugal force reduces \h\ max ^ for fib <: 10%. 
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Fig. 13. Precollapse central angular velocity Q. c ,■ versus the rota- 
tion rate rate /3b or the mass of the inner core M,- C] £ at core bounce 
for models run with the El or ST EoS. We define the mass of 
the inner core as the mass enclosed by the entropy maximum at 
core-bounce, i.e. the unshocked region. 



tonisation, drives co nvection inside the lower layers of the PNS 
dDessart et al.l l2006h at a radial position of ~ 10-30km and 
a density range of ~ 10 12 - 10 14 gcrrT 3 and therefore causes 
now the entire GW emission. The cooling PNS contracts with 
time, which causes the convective zones to migrate towards 
smaller radii and shrink. However, we point out that our leakage 
scheme overestimates the neutrino cooling processes, as shown 
in Fig. |2 Hence, this mechanism proceeds too quickly for model 
R1E1CA L (Ott 2009, private communication). The PNS convec- 
tion exhibits GW emission of roughly ~ 0.5 - 1cm amplitude, 
as can be seen in Fig. [8] (? > 20ms). The corresponding spec- 
tral distribution is shown in Fig. [10] A broad peak rises between 
~ 700 - 1200Hz and reflects the dynamical timescale tj vn of the 
violent overturn activity of a millisecond scale inside the PNS. 
The model's SNR for LIGO at lOkpc is again around unity. 
The high frequency tail of the spectrum (> 700Hz), which is 
present due to PNS convection, cannot contribute to the SNR as 
it lies below the current detector sensitivity. Our computed GW 
strains for PNS con vectio n agree roughly in amplitude with the 
ones found in lOttl d2009bl) for axisymmetric MGFLD models. 
However, the amount of released energy emitted is found to be 
about one order of magnitude higher compared to his simula- 
tions. This discrepancy is most likely due to the lower average 
frequency content of the GWs in Ott's model (~ 350Hz, see his 
Fig. 7), as dEcw/df °c v 2 . We suppose the reason for this mis- 
match to be the different radial location of the convectively un- 
stable region and thus the different related dynamical timescale 

tdyn • 

3.2. Rapidly rotating core collapse 
3.2.1. Core bounce 

Effects of the rotation rate on the GW signature 

Rapidly rotating progenitors = ^...4OTads _1 in our 

model set) undergo different core-collapse dynamics compared 
to the previously discussed non- and slowly rotating models. 
Conservation of angular momentum in combination with con- 
traction leads to a massive spin-up and hence oblate deforma- 



tion of the core. The collapse is halted either by pure stiff- 
ening of the EoS above nuclear saturation density, or, if rota- 
tion is sufficiently strong, by a combination of the centrifugal 
and the nuclear forces. The abrupt slowdown of axisymmetri- 
cally axisymmetrically-arranged and quickly rotating bulk mat- 
ter gives rise to rapid temporal variations in the quadrupole 
tensor, resulting in the emission of GWs. Note that the core 
remains essentially axisymmetric during the collapse and the 
early postbounce times ( t — tb < 10ms), as already pointed 
out in lOtt et all (l2007albh . Within the chosen parameter space 
of the rota tion rate, all our models exhibit a so-called type I 
GW burst dZwerger & Miillerl d 19971) . see Fig. Qj} around core 
bounce, no matter what the initial choice of the EoS or the 
magnetic field configura tion. This was previously also found by 
Dimmelm eier et ail d2008l) in 2D GR simulations without mag- 
netic fields. 

The question now arises what kind of information could 
possibly be delivered from a quasi-axisymmetric type I GW 
burst, since it is a priori unclear how degenerate it is with re- 
spect to the model parameters such as the EoS, rotation rate 
and so forth. This was al ready investigated in great detail by 
iDimmelmeier et al] d2008l) who performed an extensive set of 
2D GR core-collapse simulations. Our 3D results show the same 
systematics: the peak amplitude \h\ nwlx j, = A + n/R scales about 
linearly with /3b for models up to a moderately rapid rotation 
(J3b < 10%, see Fig. fT2l. The Fourier-transforms of the bounce 
wave trains (± 5ms relative to core bounce) show for most mod- 
els in the indicated parameter a spectrum with a narrow band- 
width, peaking around ~ 800 - 900Hz (see Tab. 0. Moreover, 
with growing rotation rate, prompt convective overturn in the 
rotational plane is suppressed by the influence of positive an- 
gular momentum gra dients. This effect was first pointed out by 
lEndal & Sofial (fl978l) and is known as the Solberg-H0iland in- 
stability criterion. 

The outcome of our 3 D models confirms the statement of 
Dim melmeier et ail (l2008h that two parameters are essential for 
the behaviour of the GW amplitude around t — % = 0, namely the 
mass of the inner core at bounce, which we denote as M, c .£, and 
the initial central rotation rate Q Ci ;. Over the parameter range 
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Fig. 11. Left: Time evolution of the GW amplitude A +II of the models R3E1AC (solid line), R3E2AC (dashed line) R3E3AC (dotted 
line) and R3STAC (dashed-dotted line). Right: Model R5E1 AC's GW amplitude A +II . 



covered by our models, fib is a strictly monotonic function of 
the initial central angular velocity Q c l , as displayed in Fig. [T3] 
The mass of the inner core Mfcj, is linked to fib via its depen- 
dence on Q f , (see Fig.[T3l). The positive mass offset of the inner 
core (which is approximately constant for all rotational configu- 
rations) that occurs when switching from El to ST is interpreted 
as follows. In a static initial configuration, the mass of the in- 
ner core is proportional to the square of the electron fra ction 
Y e and the entropy per baryon dGoldreich & Weberlll980h . As 
the minimum of Y e (p) for El appears at ~ 0.276 compared to 
~ 0.293 for the Shen EoS (see Fig.Q}, we attribute the mass dif- 
ference primarily to this relative difference, and secondarily to 
changes in the specific entropy which occurs as the LS EoS per- 
mits more efficient electron capture. However, note that we may 
overestimate the spread of the inner core mass M; Cj j in depen- 
dence with rotation compared to simulations which are carried 
out with full neutrino transport (Janka 2009, private communi- 
cation). The stronger the rotation becomes (at Q c> ; > 37r), the 
increased centrifugal forces start to play a dominant role, slow- 
ing down the entire dynamics of the collapse and causing the 
core to rebound at sub- or just above supra-nuclear matter den- 
sities. The imprint of such behaviour is found in the GW signa- 
ture by a smaller maximum amplitude and lower peak frequency 
compared to slowe r rotating models, as sho wn in Fig. QT| and 
Tab. [3] Similar to Dimmelm eier et alj d2008f) . we also find that 
hmax depends sensitively on the competition of both the amount 
of imposed quadrupolar deformation due to rotation, and on the 
other hand on the average density level in the for GW emission 
dynamically relevant region of the inner core. The density in 
the central region of the PNS is lowered considerably by cen- 
trifugal forces, however, no longer compensated by a prominent 
quadrupolar deformation which results from rapid rotation. The 
'optimal' configuration for strong GW emission now is overshot, 
which causes a smaller maximum amplitude and a lower peak 
frequency. 

Our findings sta n d in v ery good qualitative agreement with 
iDimmelmeier et all d2008l) . who recently performed a large 
set of 2D core-collapse simulations in GR, nearly identical 
micro-physical input, but without magnetic fields. However, 
there are some quantitative differences. For our models that 
undergo a pressure-dominated bounce, we find spectra which 
peak in average some 100-150Hz higher than the models in 



Dimmelm eier et ail (l2008h . It has been shown by Dimmelmeier 
(2007, private communication) that the difference stems from 
the fact that fully relativistic calculations shift the GW bounce 
spectrum to lower frequencies in comparison to the ones us- 
ing an effective, spherically symmetric gravitational potential. 
Furthermore, for comparable precollapse rotational configura- 
tions, our models return higher peak GW values. We suspect 
the size of the inner core is the major cause of this difference. 
The mass of the inner core for all our simulations is r oughl y 
~ O.IMq bigger than the ones of IDim melmei er et al.l ((2008). 
If we take i nto account tha t we a re using different ele ctron 
capture rates (|Bruennl d"j"985) versus Lang anke et al.l (120031) and 
iHix et al] d2003l) ). we consider the mismatch to be understood, 
as updated rates cause the mass of the inner core in our models 
to shrink. 

Effects of magnetic fields on the GW signature 

The presence of magnetic fields in our models slows down the 
accretion of angular momentum onto the PNS via field wind- 
ing. For example, the poloidal field's stress acts on fluid parti- 
cles moving in the x-y plane in a direction opposing the motion, 
leading to a deceleration. Thus, while the GW signal from the 
initially 'weakly' magnetised model R4E1 AC is already strongly 
affected by centrifugal forces at bounce (the frequency peak at 
bounce is at ~ 385Hz, while the central angular velocity is Cl c j, ~ 
6200rads~'), the initially more strongly magnetised but other- 
wise comparable model R4E1FC still undergoes a pressure dom- 
inated bounce with a frequency peak at ~860Hz and a central 
angular velocity of Q Cj j ~ 5600rads~'. However, note that this 
effect only gets prominent for initial magnetic fields that are by 
two or ders of magnitude stronger than suggested bv lHeger et al.l 
d2005l) . We will discuss the issue of strong magnetic fields in 
more detail in sec. 13.2.31 



Effects of the EoS on the GW signature 

In order to understand the dependence of the GW burst at bounce 
on the EoS, we repeated several simulations changing only the 
EoS while keeping the other parameters fixed. The most promi- 
nent change occurs when switching from El to ST. Applying the 
latter EoS in our models leads to systematically larger absolute 
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Fig. 14. Radial profiles of the weighted density pr 2 of the spher- 
ically averaged data at time t — relative to core bounce for 
model R3E1AC and its opposing model R3STAC. 



GW amplitudes at lower frequencies compared to its counter- 
parts, as shown in Tab. [3] and the left panel of Fig.Qj] This was 
also observed e.g. by iKotake et al.l d200 4b). where the two EoS 
were compared. The shift to lower frequencies can be explained 
by the fact that the typical timescale of the GW burst at bounce is 
given by the free-fall timescale Td yn ~ 1/ ^jGp ic , where p ic is the 
mean density of the inner core. Since El leads to substantially 
higher central densities at core bounce compared to ST in simu- 
lations that are not dominated by centrifugal forces, the spectral 
peak of the GW signal is shifted to higher frequencies. The max- 
imum GW burst amplitude depends on the dynamical timescale, 
the mass of the inner core as well as on the global density distri- 
bution inside th e core, as pointed out e.g. in lKotake et al .(2004b) 
and IDimmelm eier et al.l d2008l) . The dimensionless GW ampli- 
tude is roughly proportional to M, c , divided by the square of the 
dynamical timescale (h oc Mk/tJ ). Hence, it scales approxi- 
mately linearly with density. This implies that one could expect 
higher GW peak amplitudes from the more compact cores in the 
case of El. However, since models using El are slightly more 
compact in central regions, they exhibit lower densities in the 
outer l ayers. This is displayed in Fig. [7] Following the reason- 
ing of Dimmelmeier et al. (2008), we display the quantity pr 2 
in Fig. [14] It is the essential quantity in the integrand of the 
quadrupole GW formula (see Eq. ITTb . From this plot it is ap- 
parent that models run with El have higher pr 2 at small radii, 
ST yields higher values at intermediate and large radii. Due to 
their larger volume, these regions contribute more to the total 
quadrupole integral. For fast rotators (Q. c j > 3n), centrifugal 
forces sta rt to play a dominant role, as already previously ob- 
served by Dimmel meier et all d2008l) . Here, the relative differ- 
ence between the GW signatures of models run with two EoS 
decreases, since in the regime of lower densities, the EoS do not 
differ significantly. When comparing, e.g., models R3E1 AC with 
R3STAC, the absolute size of the burst amplitudes vary roughly 
25%, whereas R4E1 AC and R4STAC are only discriminated by 
~ 4% (see Tab. [3). In summary we state that it seems very dif- 
ficult to reveal information about the different two EoS by con- 
sidering the GW signature from core bounce alone. For models 
run with either the LS or the Shen EoS and all other parameters 
being indentical, the differences brought about by the EoS are 
clearly distinguishable. However, since a small variation in one 
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Fig. 15. Angular velocity profile along the positive x-axis for dif- 
ferent models at 10ms after bounce. Note that the hump in the 
angular velocity profile at about ~ 1 1km seems to be a generic 
feature, caused by accreting material with rather high specific 
angular momentum which accum ulates on t he nucl ear density 
region of the PNS, as discussed in lOtt et all d2006bl) . It is most 
pronounced in models with an initial rotation rate < R3 and gets 
flatter for high initial rotation rates. 

of the other parameters can easily have a similar effect as the 
EoS change, it will be nearly impossible to constrain the nuclear 
EoS in the general case. 

When changing the LS compressibility from K=180MeV to 
K=220MeV or K=375MeV, the features of the collapse dynam- 
ics and the corresponding GW emission remain practically unal- 
tered, as displayed in the left panel of Fig. Q~T] The only notable 
difference occurs in the vicinity of core bounce at the center of 
the PNS: Models run with the softest version of the LS EoS (El) 
allow the core to bounce at slightly higher central density com- 
pared to E2 and E3 (see Tab. [TJ. At the same time, models run 
with El exhibit lower densities at larger radii compared to cases 
where E2 or E3 was applied (see Fig. |7). However, since the 
differences in the radial density profiles are relatively low, these 
effects cancel each other if we consider again the quantity pr 2 . 
This leads to GW amplitudes of similar size and frequencies. 

3.2.2. Gravitational waves from the nonaxisymmetric 
rotational instability 

General remarks 

Rotating proto-neutron stars can be subject to non-axisymmetric 
rotational instabilities in situations when T/\W\a yn = P exceeds a 
certain critical value. Since the growing instabilities carry the 
object's spheroidal- into a triaxial configuration with a time- 
dependent quadrupole moment, strong GW emission is to be 
expected. The best understood type of instability is the clas- 
sical dynamical bar mode rotational instability with a thresh- 
old valyi£_oJ_^ HL ^^%^J^owever, strong evidence was found 
by IDimmelmeier et al.l d2008l) that it is unlikely that the PNS 
reaches rotation rates required for it to become unstable dur- 
ing the core-collapse of and the early postbounce phase the 
iron core. Another possibility is the secular instability, triggered 
at moderately high f3 sec ~ 14% if a dissipative mechanism is 
present. It grows on the relat ively slow dis sipative timescale of 
the order of a second (see, TassouJ (ll978l) V Since none of our 
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Fig. 16. Snapshots of the vorticity's z-component w z — (V x v) z in the equatorial plane for models R4STAC (upper panels, t - % — 
10, 29, 63ms) and R4E1FC/, (lower panels, t — t/, = 10, 29, 54ms) at three representative instants of their evolution. The innermost 
300 2 km 2 are displayed, and the colour is encoded in units of [s~'], ranging from -5000 (white) to 5000 (black). The upper panels of 
vorticity plots show nicely a two-armed m = 2 spiral pattern, while the middle plot of the lower panel mainly shows the same for a 
m — 1 mode. 
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Fig. 18. Spectral energy distribution of the GW signal emitted 
along the polar axis from model R4STCA at a distance of lOkpc. 



models reach such high /3-values, both instabilities cannot play 
any role in our simulations. 

However, recent work, some of which h as been car- 
ried out in idealised setups and assumptio ns ( Shibata et al. 
2002HSaiio et alJ2003tlSaiio & Yoshidal2006tlWatts et alJ2005t 
Ou & Tohlind 120061; ICerda-Duran et all 120071) and later_also 
in more self-consistent core-co llapse simulations (lOtt et al.l 
2007b|a|; IScheidegger etalJ l2008). suggests that a differentially 
rotating PNS can become dynamically unstable at 77|W|-values 
as low as ~ 1%. Today this so-called low 'T/\W\' instability is in- 



terpreted as being a resonance phenomenon dWatts et al.l l2005). 
The underlying mechanism is suspected to be the amplification 
of azimuthal (non-axisymmetric) modes at co-rotation points, 
where the pattern speed cr p = cr/m of the unstable mode matches 
the local angular velocity, (oc exp[/(<xf - m®)]) where <x is the 
mode's eigenfrequency. The PNS is differentially rotating out- 
side a radius of ~ 10km, as displayed in Fig. [15] This differen- 
tial rotation provides a reservoir of shear energy may be tapped 
by the instability. The latter leads to spiral waves, as displayed 
in Fig. [IB] which transfer angular momentum outwards (see, 
e.g., (lLovelace et al.l[l999l) ). Note that this entire phenomenon 
appears to be closely related to the Papaloizou-Pringle instabil- 
ity w hich occurs in accretion disc s around a central gravitating 
body (fPapaloiz ou & Pringlelll985l) . 

In order to investigate the growth of the non-axisymmetric 
structures, we monitor the PNS by decomposing the density at a 
given radius R in the equatorial plane (z = 0) into its azimuthal 
Fourier components: 



p(R,Z,cf>)= ]T C m (R,z)e im * 



1 C 2n 

In Jo 



p(R,z, 



(31) 



(32) 



The gravitational wave morphology resulting from the non- 
axisymmetric process generally shows narrow-band and highly 
periodic signals which persist until the end of our simulations. 
This is shown in Fig.fTTIand the upper panels of Fig.|2T| Bearing 
in mind that the effectively measured GW amplitude scales with 
the number of GW cycles N as h eff oc h V^V (see lThornel (fl989)) 
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Fig. 17. The upper left panel shows the time evolution of the GW amplitudes A + and A x of model R3E1 AC emitted along the polar 
axis. The upper right displays the same for model R3STAC. The lower two panels contain the same information for R4E1 AC and 
R4STAC from left to right. Note that A + and A x oscillate at the same frequency, phase shifted by n/2. 



and that it could last for several hundreds of ms as our simu- 
lations suggest, the chances of being able to detect such kind 
of signal are enhanced compared e.g. to the short duration burst 
signal at bounce, as one can see in Fig. [18] In the lower pan- 
els of Fig. [2J] the normalised mode amplitudes of the models 
R3E1AC L and R4E1FC L , with 

A m = \C m \/C = -L V^(C„,) 2 + 3(C,„) 2 , (33) 

are plotted in order to follow the behaviour of unstable modes. 
We generally find modes with density wave numbers m = 1, 2, 3 
being triggered, with the m = 1 or m — 2 as the overall dominant 
ones, depending on the individual model. Furthermore, note that 
all modes have the sam e pattern speed, as previou sly observed by 
lOtt et al.l (l2007blah and lScheidegger et al) d2008l) . In Fig. [16] the 
upper panels of vorticity plots show nicely a two-armed m — 2 
spiral pattern, while the middle plot of the lower panel mainly 
shows the same for am- 1 mode, as we expect from a mode 
analysis of the corresponding model R4E1FCl presented in Fig. 
|2T1 After the early linear growth phase, the modes saturate due 
to Kelvin-Helmholtz shear instabilities, which break the spirals 
apart in the outer layers, as dis played in Fig.fTBIand prev iously 
observed and discussed in e.g. Cerda-Duran et al. (2007). Note 
the close relation between the m = 2 bar mode and the emis- 
sion of GWs. The growth and saturation of this mode is im- 
printed on the GWs emitted. As soon as it exceeds the Cartesian 
m — 4 noise background, strong GW emission at a frequency 



corresponding to twice the m — 2 pattern speed along the pole 
emerges, with the + and x polarisations shifted by a quarter cy- 
cle, as one could expect from GWs emitted by a spinning bar. 

Models without deleptonisation in the postbounce phase: 
Effects of the EoS and magnetic fields on the GW signature 

We now turn the discussion in more detail to individual sim- 
ulations. Models with an initial rotation rate of at least R3 
(fib ~ 5.2%) sooner or later become low (5 dynamically un- 
stable in our parameter set (see Table O. Note, however, that 
IScheidegger et al.l (|2JX)9j) recently found models to become dy- 
namically unstable at even slower initial rotation rates (fib ~ 
2.3%). 

When comparing models that were carried out with the 
Shen- and the LS (El) EoS, we find that the 'Shen'-models emit 
GWs at significantly higher frequencies than their LS counter- 
parts. This result is a consequence of the fact that the specific 
angular momentum, whi ch scales as j 7 <x pu> z , is roughly pre- 
served on a mass shell (Ke il et al.lll996l) : The innermost part 
of the PNS in 'Shen' -simulations, which rotates nearly in per- 
fect solid body rotation at the pattern speed (see Fig. [19), has 
lower density in this radial region. Hence the central rotation 
rate must in turn be higher to fulfill the conservation law. As 
an example, we compare the models R3E1AC and R3STAC in 
detail. While R3E1AC rotates at t - % = 10ms with a cen- 
tral rotation rate of lorjeiac ~ 2300rads~ 1 and has a cen- 
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Fig. 19. The upper panel shows the angular velocity profile at t - 
tb ~ 10 ms along the positive x-axis for the models R3E1 AC and 
R3STAC. The pattern speeds of the corresponding simulations 
are indicated by horizontal lines. The lower panel displays the 
spherically averaged radial density profiles at the same particular 
time. 
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Fig. 20. The figure shows the enclosed mass [M Q ] as a function 
of radius for models R4E1FC and R4E1FC L at t - t b -25 ms, 
which explains the strong GW emission of the latter model. 



PNS. The detailed discussion of this issue is postponed to sec. 
137231 

Centrifugal forces set a limit to the maximum frequency of 
the GW signal in similar fashion to the situation at core bounce. 
As discussed in the previous subsection, the limit is somewhere 
around ~ 935Hz (which is twice the pattern speed!). The faster 
the initial rotation rate, the stronger the influence of centrifugal 
forces, which slow down in the postbounce phase the advection 
of angular momentum onto the PNS. The result is a slower cen- 
tral rotation rate, a lower pattern speed and thus GW emission 
at lower frequencies. Beside these semi-quantitative statements 
which allow for the distinction of the simulations' input physics 
by its GW signature on a model-to-model basis, where only one 
parameter is varied while keeping the others fixed, it is in general 
very difficult to discern effects of individual features of the input 
physics in a GW signal that cannot unambiguously be attributed 
to a specific model. The degeneracy in the simulation results is 
large with respect to the rotation rate, the magnetic fields and the 
underlying EoS. 



tral density p RiS TAC ~ 3.38xl0 I4 gc m ~ 3 > R3STAC revolves with 
the values of ojristac ~ 2800 rads~' and Prieiac ~ 2.63 x 
10 14 gcm~ 3 , as displayed in Fig. |T9j Doing the maths, the ratios 
of cjrieiac/lor3stac top R 3STAc/PR3STAC are about the same. We 
also state that the dynamical instability in the 'Shen' cases grows 
faster generically than in the LS simulations, as shown in Fig. 
El 

We generally observe slower growth of the T/|W| unstable 
modes in situations where we applied stronger initial poloidal- 
than toroidal magnetic fields. Although such magnetic fields 
(Bpoij > B tor .i) may not be motivated by stellar evolution cal- 
culations (iHeger et al.ll2005h . it is still important to study their 
effects. Model R3E1AC for example starts to emit strong GWs 
due to the low T/\W\ instability around ~ 50ms after bounce, 
while model R3E1DB does not within the duration of the simu- 
lation, which we followed until ~ 65ms after bounce (note, how- 
ever, that the latter model shows strong growth of the m=l,2,3 
modes although GW emission due to the low T/\W\ did not set 
in yet). The poloidal fields are able to suppress the dynamical 
instability for some time as they slow down the spin-up of the 



Models with deleptonisation during the postbounce phase 

The improved input physics of the two leakage models 
R3E1ACl and R4E1FCl allows us to address the question how 
the inclusion of deleptonisation in the postbounce phase quan- 
titatively alters the GW signal of our 3D MHD models. While 
the results show qualitative agreement with our earlier findings, 
they clearly deviate in quantitative terms. As the most strik- 
ing feature, these models show 5 - lOx bigger maximum GW 
amplitudes due to the nonaxisymmetric dynamics compared to 
their counterparts that neglect neutrino cooling, as one can see 
in Table [3] or when comparing Fig. [T71 with Fig. [21] This sug- 
gests that the treatment of postbounce neutrino cooling plays an 
important role when it comes to the quantitative forecast of GW 
signals from a low B instability. The neutrino cooling during the 
postbounce phase leads to a more condensed PNS with a shorter 
dynamical timescale compared to the purely hydrodynamical 
treatment, as shown in Fig. [20] This in turn is directly reflected in 
the dynamical evolution: The shock wave stalls at considerably 
smaller radii and becomes more quickly unstable to azimuthal 
fluid modes (see Fig. [TBI ). Since there is much more matter in 
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Fig. 21. The upper left panel displays the emission of the A + and the A x -amplitude along the pole for model R3E1AC/.. The upper 
right panel displays the same for model R4E1FC/,. The lower two panels show for the same models as above the normalised mode 
amplitudes A,„ for m = 1, 2, 3, 4 extracted at a radius of 25km. 



the unstable region of these models, the unstable modes grow 
faster and the triggered spiral density waves cause the emission 
of much more powerful GW. 



We close this subsection by pointing out that p revi- 
ous core-collapse comp utations by lOtt et alJ d2007blah and 
Scheidegger et al. (2008) show both qualitative and quantitative 
agreement with the ones presented in the previous subsection as 
they incorporate nearly identical micro-physics. However, they 
mismatch on a quantitative scale with R3E1ACl & R4E1FCl 
due to the absence of the postbounce neutrino treatment. 
However, we point out that our leakage scheme overestimates 
the compactification of the PNS due to neutrino cooling. The 
'reality' for the strength of GW emission therefore sould lay in 
between the results from the pure MHD and the leakage treat- 
ment. We also want to point out that another limitation which 
might affect the absolute values of the GW signature from a 
T/\W\ dynamical instability is the grid resolution. Our choice 
of a uniform grid leads to more than sufficient resolution at the 
stalled sh ock, but may under r esolv e the surface of the PNS. For 
example. ICerda-Duran et all (120071) showed that resolution can 
have a significant effect on the instability's developement and its 
GW signal. 



Ertropy [kB/baryon]. time-tb = 0.01 8B48 




Fig. 23. Snapshot of model R4ElCF's entropy distribution in the 
first octant at a representative instant of its evolution. The inner- 
most 300 3 km 3 are displayed. 
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Fig. 22. The upper left panel shows the first 20ms of the time evolution of the rotational energy parameter B = T /\W\ of models 
R4E1CF, R4E1FC and R4E1FCl- The upper right panel displays the time evolution of the magnetic energy parameter B mag for the 
same models. The small hump in the f3„ mg profiles at t-tb = 0, best visible for R4E1FC and R4E1FC/,, stems from the fact that field 
compression scales as B oc p 2 ^ . When the core overshoots its equilibrium position at bounce and re-expands, the density decreases 
temporarily, which in turn is reflected in the energy density of the magnetic field. The lower left panel shows the evolution of the 
magnetic parameter B lormag , B po ^„ mg and the total B mag for R4E1CF, the lower right displays the same quantities for R4E1FC. 
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Fig. 24. Model R4ElCF's time evolution of the total quadrupole 
amplitude A+n,tot (left panel) and its contributions from the mag- 
netic field A + n, mag (right panel). 



3.2.3. The influence of strong magnetic fields on the 
gravitational wave signature 

In the previous three subsections we have shown that the core- 
collapse dynamics and thus the GW signal of initially weakly or 
moderately strong magnetised stellar cores (< 10 n G) is hardly 
affected by magnetic fields. However, in the case of strong initial 
fields (~ 10 12 G) things change dra matically. The com bined ac- 
tion of flux-freezing, field winding dMeier et al.lll976l) and also 
of the magneto-rotational instability MRI ( Balbus & Hawlevl 



1998) generally may lead to growth of magnetic fields by 
many orders of magnitude to values where the magnetic pres- 
sure reaches the order of magnitude of matter pressure. This, 
in turn, trigge r s a co llimated, bipolar jet-explosion (see e.g. 
iBurrows et al.l d2007l) . and references therein) by converting 
magnetic energy into kinetic energy. However, with our current 
grid setup we are unable to resolv e exponential growth triggered 
by the MRI dEtienne et alj|200q) ; the initial magnetic fields are 
solely amplified by compression during the infall phase and by 
magnetic winding. 

For two reasons, magnetically-driven explosions are of in- 
terest for the prediction of GWs. Fir stly, the bipola r outflow of 
matter results in a 'memory e ffect' (jThornd 19891) in the GW 
signature, as observ ed e.g. by [Obergaulingeret al- (200$ and 
Shib ata et alJ (120061) . As the out-stream of matter usually hap- 
pens along the rotational axis, the amplitude A + n will grow over 
time, scaling as A+n K 2;«v?, where the ejected mass m increases 
constantly in the early stage of a magneto-rotational supernova. 
Secondly, the GW amplitude of these (not necessarily) realistic 
models is also affected passively by strong magnetic fields, as 
they rise during collapse and early postbounce phase to values as 
high as ~ 10 16 ~ 17 G. The magnetic energy density may provide 
sizable contributions to the overall amplitude, as estimated ap- 
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proximately by the following formula (cf. iKotake et al.1 d2~004b) 
Eq. 22): 



Tjmag 



= ~^-r ~ 10% ( 

PcC 1 \ 



several x 10 



17 g) 



10°gcm- 



(34) 



where B c and p c stand for the central values of the magnetic 
field and density. If magnetic contributions become significant it 
is necessary to include them in Eq.[23] which then turns into Eq. 

m 

For the purpose of studying the effect of very strong mag- 
netic fields on the GW signal in 3D, we have carried out three 
runs: R4E1CF, R4E1FC and a leakage model R4E1FC L . In order 
to overcome the technically chall enging simulation of mag netic 
field growth through small scale dObergaulinger et aLll2009l) or 
long-term processes, we applied initial configurations which 
may not be realised in nature, but can deliver the formation of 
jets in a similar way as natural field growth would. 

Model R4E1CF, which was set up with very strong initial 
poloidal fields (see Table [TJ, shows a typical behaviour for a 
magneto-rotational core collapse followed by a jet-like explo- 
sion. Field compression during the collapse phase strongly am- 
plifies both the toroidal- and poloidal magnetic fields, since flux- 
freezing in stellar collapse guarantees the B-field to scale as 
B ~ p 1 ^ . Note that we evaluate toroidal and poloidal magnetic 
field components as: 



B„ 



-V* 



+ Bl 



B 



pol 



B- 



(35) 
(36) 



Furthermore, the toroidal field is build up by tapping the mag- 
netic energy contained in the poloidal component of the field 
trough winding, whilst the poloidal field component grows only 
little during the postbounce phase by action of meridional mo- 
tions in the core as shown in the lower left panel of Fig. [22] In 
turn, the generated hoop stresses grow fast to reach high values 
near the polar region, and the ratio of magnetic to matter to fluid 
pressure reaches P magi P matter ~ 1. A bipolar collimated jet ex- 
plosion is launched, quickly leaving the computational domain 
(see Fig. 1231. 

Another remarkable feature is that, in contrast to weakly 
magnetised models, magnetic breaking efficiently decelerates 
the in ner core by redistributing angular momentum (Meie r et alj 
Il976l) . This is shown in the upper left panel of Fig. [221 Note that 
model R4E1FC slows down at a higher rate than R4ElFCz. as 
its fluid can re-expand to larger radii due to the absence of neu- 
trino cooling. The GW signal is displayed in Fig. [24] In order to 
distinguish the different contributions to the total GW amplitude 
hJj T tol , we split the magnetic and fluid part in the following way: 



h" =h" +h" 

,l ij,tot ' l ij,matter T ' l ij,nuig- 



(37) 



Around core bounce, the structure of the GW signal is very sim- 
ilar to that of less magnetised cores, exhibiting a clear type I 
signature. However, shortly after bounce we observe a growing, 
slowly time-varying offset of A + n relative to the horizontal axis 
compared to ring-down oscillations of the weakly-magnetised 
simulations. There are two reasons for this behaviour. Firstly, the 
magnetic contribution A + n, m ag to the total GW amplitude grows 
strongly as magnetic forces act on the core. Secondly, as we 
have explained earlier in this subsection, matter-outflow along 
the z-axis also contribut es to the signal . This b ehaviour was al- 
ready observed by lObergaulinger et all d2006l) in axisymmetric 
simulations and its characteristical GW amplitude was named 



'type IV signal' (Obergaulingeret al. 2006). The contribution of 
the magnetic amplitudes shrinks with the onset of the jet. The 
emerging matter gains its kinetic energy by tapping the energy 
stored in the magnetic field, which causes a drop in B mag and 
hence also in A+n, m ag, as it can be seen in Figs.[22]and[24]at about 
t — tb ~ 10ms. The non-axisymmetric amplitudes are negligible 
compared to the axisymmetric part of the wave train (see Tab. 
[3} and the result of some prompt convective motions and com- 
plicated time-variations in the postbounce magnetic field con- 
figuration. However, even if a Galactic supernova was optimally 
orientated for the detection of a type IV signal, it is still ques- 
tionable whether we could distinguish it in early stages from an 
ordinary type I signal. The charact eristic offset of this particu lar 
signal type was suggested e.g. by Obergaulinger et al.l d2006l) to 
be a valid measure of aspherity of the ongoing supernova explo- 
sion. However, since the memory effect in the amplitude appears 
on the long timescale of several times 10ms, it would be out of 
the LIGO band. However, the pla nned space-based DECIGO in- 
strument (Kawamu ra et alJ l2006) could permit to track the low 
frequency contribution of such a GW signal in the future. 

We point out that at the onset of jet formation around t-% w 
10ms, the absolute value of the magnetic field in the polar region 
at the edge of the PNS is ~ 10 I6 G, which translates into a mildly 
relativistic fast Alfven speed. Additionally, the velocity of the 
ejected matter accelerates up to radial velocities of ~ 0.1c when 
leaving the boundary of our computational domain. These two 
points challenge the Newtonian treatment of the dynamics in our 
scheme and suggest to use a special or general relativistic code 
for the further evolution of the jet dynamics. 

Models R4E1FC and R4E1FC L show a very different dy- 
namical outcome. For these two simulations, we assumed the ini- 
tial toroidal field to be 10 3 times the values of the poloidal field, 
as sug gested from stellar evolution calculations by iHeger etail 
d2005l) . During collapse the magnetic field components primar- 
ily grow due to compression as one can see in the lower right 
panel of Fig. [22] Right after bounce, the strong toroidal fields 
in both models cause the onset of a jet, as one can see in the 
left panel of Fig. [25] However, surprisingly and in c ontrast to the 
previo usly discussed models, or simulations from Kotak e et al.l 
(2004b) where they applied similar initial conditions and ob- 
tained jet explosions in 2D, the wind-up of the poloidal- into 
the toroidal field does not occur efficiently enough. The spiral 
wave of the standing accretion shock instability (SASI, see e.g. 
iBlondin et al.l(l2003lV ). which forms at the same time, hinders and 
delays the growth of a jet, as displayed in the right panel of Fig. 
|251 We interpret this phenomenon as follows: Matter can easily 
slip along, but not move perpendicular to magnetic field. In the 
x-y plane, the poloidal field's stress acts on matter in the opposite 
direction of the fluid motion. This is also reflected in the upper 
left panel of Fig. [22] It displays clearly that the magnetic forces 
are not capable of slowing down the inner core in case of strong 
initial toroidal fields as effectively as the poloidal ones that are 
anchored in the outer stellar layers. Since models R4E1FC and 
R4E1FCl have relatively weak initial poloidal magnetic fields, 
the fluid, which rotates around the z-axis, can develop nearly 
unhindered instabilities. The developing spiral waves then coun- 
teract successfully the formation of a jet-like explosion by turn- 
ing matter aside the pole (see Fig.EBTl. where the magnetic hoop 
stresses are strongest (and in axisymmetry most probably would 
be able to launch a jet at this stage of the simulation). This cor- 
responds also to our observation from the last subsection where 
we stated that the low instability grows slower in the pres- 

ence of dominant poloidal fields, because they cause stresses that 
act against the spiral instabilities. However, further discussion of 
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Fig. 25. Snapshots of model R4E1FCl's entropy distribution in the first octant at two representative instants of its evolution. The 
innermost 300 3 km 3 are displayed. 



this phenomenon is beyond the scope of this work and will be in- 
vestigated in a subsequent study. 

The resulting GW signals consequently show a type I sig- 
nal at core bounce, subsequently followed by a low (3 instability 
(see the right panels of[2Tb. Although the GW contributions due 
to magnetic stresses in model R4E1FC/. show qualitatively the 
same features as the ones from model R4E1CF, they are smaller 
and dominated by the hydrodynamics part of the amplitude (e.g., 
^4-+ii,mag ~ 10cm vs. A+n,matter ~150cm at 30ms after bounce). 

4. Summary and Conclusions 

Core-collapse supernovae are a source of GWs that can only 
be modelled imperfectly due to large uncertainties in the ini- 
tial conditions, input physics and the technically challenging 3D 
neutrino transport. The parameter space of possible initial condi- 
tions is huge since many progenitor configurations are possible 
at the onset of collapse. In this paper we tried to outline some 
dependences of the 3D GW form upon a variety of these con- 
ditions, since for un-modelled burst analysis techiques, partial 
information on wave forms is already useful, e.g. typical wave- 
form features, approximate spectra and the type of polarisation 
and burst duration. 

With our model series containing 25 three-dimensional 
MHD core-collapse supernova simulations, we have tried to 
probe the GW signature with respect to different nuclear equa- 
tions of state, rotation rates, poloidal and toroidal magnetic 
fields, and a postbounce deleptonisation s cheme 



Similar to the findi n gs of e.g. iMuller et all J2004); 



Dimmelm eier et alJ (l2008l) : TMarek et all (l2009t) : idtil (I2009bh : 
Scheidegger et al.l ([2009), our results show that all non- and 
slowly rotating models release GWs due to prompt convection 
within the first ~ 30ms postbounce that is caused by the presence 
of a negative radial entropy gradient. Furthermore, we could 
show in simulations without deleptonisation in the postbounce 
phase that the waveforms obtained from this early stage of a su- 
pernova explosion contain indirect information about the under- 
lying EoS. Due to different radial locations of the convectively 



unstable region and the amount of matter it contains, we were 
able to distinguish the LS EoS from the one of Shen. While the 
LS EoS leads to GW emission in a frequency band peaking be- 
tween ~ 150-500Hz, the spectrum of the models using the Shen 
EoS is restricted to roughly ~ 150 - 350Hz. However, LIGO's 
current sensitivity makes it impossible to see the high-frequency 
component of the LS models. Thus, a distinction between the 
two EoS is currently not possible. Nevertheless, planned up- 
grades of the interferometers in the near future should enable 
the discrimination between the prompt convection GW signal of 
the LS- and the Shen EoS. We also found minor deviations in the 
GW characteristics for simulations which were carried out with 
different compressibility versions of the LS EoS. However, the 
differences in the frequency domain of the GW signal are negli- 
gibly small and thus not likely to be constrained by observation. 

With the inclusion of a neutrino leakage scheme in a slowly 
rotating model for the postbounce phase, we could show that the 
GWs emitted during the first ~ 20ms after bounce are predomi- 
nantly due to entropy driven 'prompt' convection. In purely hy- 
drodynamical models, the GW emission ends with the decaying 
negative entropy gradient. In more realistic models, a negative 
radial lepton gradient, caused by the neutronisation burst and 
the subsequent deleptonisation, takes over as driving force of the 
convective activity at the edge of the PNS. The long-lasting GW 
emission associated with this dynamical feature is roughly of the 
size of ~ 1cm and has a spectral peak around ~ 700 - 1200Hz. 

In our set of models, simulations with a precollapse core an- 
gular velocity within the parameter range of Q e( - = ji . . .4OTads _I 
undergo a rotational core collapse. The models all exhibit a so- 
called type I GW burst at core bounce. As the most important 
outco me, our 3D MHD models c ould confirm the recent find- 
ings of iDimmelmeier et al.l (|2008) that the purely axisymmetric 
(1=2, m=0) peak amplitude scales about linearly with the rota- 



tion rate fit, at core bounce (\h„ 



p b ) for fi b < 10%, while 



the Fourier-transform of the bounce wave trains for most models 
in the indicated parameter range align around a spectral peak of 
~ 800 - 900Hz. However, for very fast initial rotation rates of 
Q. c j > 3nmds~ l , centrifugal forces significantly decelerate col- 
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lapse and core bounce. The longer timescale and the weakened 
spin-up of the core due to the action of centrifugal forces leads 
generically to a decrease of the peak amplitude and a broad- 
ened spectral peak at lower frequencies. Furthermore, our re- 
sults indicate that the particular choice of the nuclear EoS has 
little influence on the GW signal from rotational core bounce. 
These findings are in good qualitative agreement with the ones 
iDimmelmeier et all d2008l) derived from axisymmetric models. 

Models with a rotation rate of fib ^ 5% at core bounce be- 
come subject to a so-called low T/\W\ instability of dominant 
m — 1 or m — 2 character within the first several tens of ms 
after bounce, depending on the individual model. This nonax- 
isymmetric dynamical shear instability leads to prolonged nar- 
row band GW emission at a frequency of twice the rotation rate 
of the innermost part of the PNS that rotates as a solid body. The 
fact that the effectively measured GW amplitude scales with the 
number of GW cycles N as h e /f oc h yN suggests that the detec- 
tion of such a signal is tremendously enhanced compared to e.g. 
the short-lived GW bursts from core bounce, and would allow 
us to probe the rotational state of the PNS over a long period. 
However, we point out that such a mechanism only operates if 
the progenitor is rotating much faster than predicted for most 
stars by stellar evolution calculations dHeger et al.l 12005). We 
also find that centrifugal forces set a limit to the maximum fre- 
quency of this periodic GW signal somewhere around ~ 935Hz 
as they suppress the inward advection of angular momentum. 
Besides that, we point out that GWs from a low T/\W\ instabil- 
ity are highly degenerate with respect to initial rotation rate, EoS 
and magnetic fields. Thus, it is very difficult to extract individual 
features of the input physics from the GW signal that can clearly 
be attributed to the initial conditions of a progenitor. Rapidly ro- 
tating models that include the postbounce neutrino physics at a 
qualitative level reproduce the previous findings. However, the 
GW signature from these more advanced models show huge 
quantitative deviations from the ones that treat the postbounce 
phase purely hydrodynamically. As the neutrino cooling during 
the postbounce phase leads to a more condensed PNS, the un- 
stable regions contain considerably more mass, which then re- 
sults in 5 to 10 times bigger GW amplitudes. Motivated by these 
findings, we will continue to improve the postbounce neutrino 
physics in future simulations by including neutrino heating and 
the emission of yu/r neutrinos. 

Our simulations show that the impact of magnetic fields on 
the overall supernova dynamics is generally small in cores with 
relatively weak precollapse fields (B < 10 n G). Nevertheless, if 
we impose very strong and probably unrealistic initial poloidal 
magnetic fields (~ 10 12 G), the combined action of flux-freezing 
and field winding allows the toroidal field component to grow 
over many orders of magnitude to values where the magnetic 
pressure can trigger a jet-like explosion along the poles. The 
bipolar outflow of matter then causes a type IV gravitational sig- 
nal. However, if we assume a strong initial toroidal magnetic 
field, the onset of a jet is effectively suppressed by the fast grow- 
ing spiral waves of the S ASI. This finding stands in contradiction 
to 2D simulations, where similar config urations led to 'jet'-like 
explosions, as e.g. in Kot ake et al ](|2QP4b). Thus, in simulations 
where the toroidal component of the magnetic field dominantes 
over the poloidal one, the magnetic contributions to the GW sig- 
nal are dominated by the hydrodynamical part of the amplitude. 
This effect is even stronger in simulations with a deleptonisation 
scheme in the postbounce phase. 

The major limitation of our code now is in the monopole 
treatment of gravity, since it cannot account for spiral structures, 
which could be reflected in GW. We are currently working on the 



improvement of this issue. The IDSA includes at present only 
the domi nant reactions relevant to the neutrino transport prob- 
lem (see lLiebendorfer et al.l (12009) for details). Future upgrades 
will also include contributions from electron-neutrino scattering, 
which are indispensable during the collapse phase. The inclu- 
sion of this reaction will also make the cumbersome switch of 
the neutrino parametrisation scheme to the IDSA at bounce ob- 
solete. Finally, we work on the inclusion of fi and t neutrinos, 
which are very important for the cooling of the PNS to its final 
stage as neutron star. 
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